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PREFACE 


This  report  was  prepared  by  the  University  of  Dayton  Research  Institute 
for  the  Federal  Aviation  Administration  Technical  Center  under  Contract 
FA74WA-3532  during  the  period  April  1980  through  March  1981.  The  report 
describes  the  third  refined  version  of  the  Dayton  Aircraft  Cabin  Fire  Model 
(DACFIR). 

Work  was  performed  at  the  University  of  Dayton  under  the  supervision  of 
Nicholas  Engler.  Others  at  the  Unlversi-ty  who  contributed  to  this  program  are 
John  Myers,  Steven  Vondrell,  Thomas  Shirley,  and  Zalfa  Challita.  Much  of  the 
original  development  of  the  DACFIR  model  is  due  to  Jerry  Reeves,  Peter  Kahut, 
and  James  Luers.  The  author  wishes  to  thank  Gretchen  Walther,  Jacquelin 
Aldrich,  Nahla  Abdelnour,  and  Pamela  Ecker  for  their  assistance  in  preparing 
this  report. 
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SECTION  1 
INTRODUCTION 


1.1  PURPOSE 

The  report  describes  Version  3  of  the  Dayton  Aircraft  Cabin  Fire  Model 
(DACFIR3),  a  mathematical  model  of  fire  within  the  cabin  of  a  commercial 
transport  aircraft.  DACFIR3  is  the  most  recent  refinement  of  models  developed 
to  predict  the  accumulation  of  smoke,  heat,  and  toxic  gases  in  a  cabin  by 
using  laboratory  test  data  on  the  cabin  interior  materials.  The  refined  model 
takes  into  account  compartmentation  of  the  cabin  and  the  entry  of  the  flames 
and  gases  of  an  external  fire.  In  addition,  DACFIR3  incorporates  a  number  of 
recent  advances  in  theory  and  practice  of  modeling  enclosed  fires. 


1.2  BACKGROUND 

The  DACFIR  model  was  developed  to  better  understand  the  relationship  of 
small-scale  fire  test  data  on  individual  cabin  materials  to  the  behavior  of 
these  materials  when  involved  in  an  actual  full-scale  fire.  For  many  reasons, 
this  relationship  is  not  always  straightforward.  This  fact  has  obvious 
bearing  on  fire  safety  standards  for  cabin  materials.  Mathematical  modeling 
and  computer  simulation  of  cabin  fire  development  provides  one  means  of 
clarifying  the  relationship  by  making  it  possible  to  account  for  the  complex 
geometry  of  the  cabin  and  the  various  conditions  of  ignition,  fire  exposure, 
and  ventilation. 

The  first  version  of  DACFIR,  described  in  Reference  1,  predicted  fire 
spread  and  smoke,  heat,  and  gas  accumulation  for  complex  cabin  geometries  but 
was  limited  to  wide-body  cabin  configurations.  In  order  to  compare  the  pre¬ 
diction  of  It-  -nodel  to  existing  fire  test  data  on  standard  width  cabins,  a 
number  of  refinements  and  extensions  were  made  (Reference  2)  which  resulted  in 
the  second  version,  DACFIR2.  This  second  version  accomodates  both  standard 
and  wide-body  cabins  and  predicts  oxygen  concentration  in  the  cabin  atmosphere 
In  addition  to  the  original  quantities.  Other  refinements  Included  adding  a 
capability  for  forced  ventilation  In  Improvements  to  the  methods  of  computing 
flame  and  gas  radiation. 

Recent  test  and  experiment  work  (References  3,  4,  5,  and  6)  has  focused 
attention  upon  a  highly  probable  mode  of  Ignition  for  cabin  fires— the  entry 
of  flames,  radiation,  and  gases  of  an  exterior  fire  through  a  cabin  door  or 
other  opening.  The  first  two  versions  of  the  model  addressed  only  fires  ori¬ 
ginating  within  the  cabin,  and  so  cannot  be  used  to  study  this  interesting 
case.  Also,  division  or  compartmentation  of  the  cabin  to  contain  a  growing 
fire  has  been  proposed  as  a  fire  safety  measure.  This  concept  has  been  exa¬ 
mined  experimentally  (Reference  7)  but  no  theoretical  methods  for  its  study 
have  been  available.  Finally,  a  number  of  advances  (e.g..  References  8  and  9) 
in  both  the  theory  and  practice  of  modeling  enclosed  fires  have  appeared  since 
the  original  conception  of  DACFIR  that  can  be  Incorporated  in  the  model  to 
bring  It  up  to  current  standards. 

These  three  factors— Interest  In  exterior  fire  ignition.  In 
compartmentation,  and  the  current  Improvements  in  mathematical  fire  modeling— 
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have  resulted  in  a  third  version  of  DACFIR  which  is  the  subject  of  this 
report.  The  report  Is  divided  into  two  volumes.  The  first  volume  presents 
the  assumptions,  theory,  and  formulation  of  DACFIR3  including  a  brief  descrip¬ 
tion  of  the  computer  code.  Volume  2  contains  several  appendices  on  deriva¬ 
tions  of  some  of  the  model  equations,  a  user's  guide  for  the  DACFIRS  program, 
and  a  listing  of  the  computer  code. 
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SECTION  2 

STRUCTURE  OF  THE  MODEL 


The  formulation,  assumptions,  and  structure  of  the  DACFIR  model  reflect 
the  model's  basic  purpose:  to  predict  the  accumulation  of  smoke,  heat,  and 
toxic  gases  during  the  emergency  evacuation  period.  Our  interest  is  in  the 
conditions  which  may  impede  evacuation  and  the  role  the  cabin  materials  have 
in  the  development  of  these  conditions.  Therefore  the  model  deals  with  the 
early  stages  of  the  fire;  that  is,  times  when  only  a  small  percentage  of  all 
the  cabin  materials  are  Involved.  During  this  period  the  cabin  Is  filling 
with  combustion  products  but  evacuation  may  still  be  possible.  The  rate  at 
which  the  cabin  becomes  uninhabitable  is  the  key  factor  that  the  model  seeks 
to  predict.  This  time  scale— perhaps  the  first  five  or  ten  minutes  after  a 
fire  develops  in  the  cabin— has  the  major  Influence  on  the  model's  structure 
in  terms  of  the  physical  processes  that  must  be  considered  and  how  the  pro¬ 
cesses  are  handled  mathematically. 

Because  the  model  must  attempt  to  deal  comprehensively  with  all  the 
Important  phenomena  of  the  problem,  the  level  of  detail  and  sophistication  of 
the  component  pieces  of  the  model  Is  necessarily  somewhat  simple.  However, 
when  the  simple  models  of  each  phenomenon  are  assembled  into  the  complete 
system,  the  resulting  structure  is  complex.  This  section  gives  an  overview  of 
the  complete  model  by  briefly  explaining  the  modeling  techniques  and  assump¬ 
tion  and  the  relationship  among  them.  More  detailed  explanations  of  the 
separate  components  are  given  in  Sections  3,  4,  and  5. 


2.1  CABIN  GEOMETRY  AND  FIRE  GROWTH 

The  DACFIR  model  considers  fires  originating  on  and  propagating  over  the 
fixed  interior  surface  of  the  cabin.  The  model  makes  detailed  predictions  of 
the  areas  of  each  type  of  material  on  fire  at  any  time.  The  rate  at  which  a 
fire  develops  depends  upon  the  types  of  materials,  the  orientation  of  the 
burning  surface,  and  the  thermal  conditions  in  the  cabin.  We  approximate  the 
cabin  geometry  as  being  constructed  of  only  horizontal  or  vertical  surfaces. 
This  assumption  simplifies  the  calculations  for  simulating  the  fire  growth 
while  retaining  the  Important  geometric  features.  The  method  of  predicting 
the  burning  area  on  each  surface  uses  a  gridding  scheme  which  divides  the  fuel 
surface  Into  discrete  regions.  Each  surface  within  the  cabin  Is  divided  Into 
equal  size  squares  using  a  grid  of  vertical  and  horizontal  lines.  A  con¬ 
venient  dimension  for  the  square  size  Is  six  inches  (0.5  ft)  on  a  side.  These 
fuel  squares  are  referred  to  as  "elements"  -In  the  model.  Figure  1  Illustrates 
how  the  cabin  Interior  Is  approximated  In  this  method.  Seclton  3  presents  the 
details  of  the  gridding  scheme. 
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combustion  of  solids  does  not  let  us  obtain  the  pertinent  flammability  charac¬ 
teristics  from  first  principles  for  the  complex  materials  and  structures  of  an 
aircraft  cabin.  However,  these  characteristics  may  be  measured  directly  in 
laboratory  experiments.  The  DACFIR  model  has  been  developed  assuming  that 
these  characteristics  can  be  supplied  as  input  data.  The  specific  material 
flammability  information  used  by  DACFIR  is  shown  in  Table  1. 

These  flammability  characteristics  are  not  intrinsic  material  properties 
but  are  also  functions  of  exterior  conditions  during  burning.  For  DACFIR  we 
have  made  the  simplifying  assumption  that  these  items  are  functions  only  of 
the  externally  applied  radiation  flux.  While  this  is  certainly  an  over¬ 
simplification,  since  other  factors  such  as  fraction  burned  or  sample  orien¬ 
tation  are  definitely  important,  current  opinion  (e.g.  Reference  10)  holds 
that  flame  radiation  has  the  dominant  effect  on  burning  rate  for  fires  of  the 
scale  considered  by  DACFIR3. 

Table  1  contains  two  emission  rates  for  each  material,  one  for 
smoldering  and  the  other  for  burning  (open  flaming).  The  smoldering  condition 
is  initiated  when  the  heat  flux  to  the  element  reaches  a  specified  level  and 
the  material  decomposes  without  flames.  Smoke  and  toxic  gases  are  emitted  at 
a  constant  rate  until  (1)  the  heat  flux  is  reduced  below  the  specified  level, 
(2)  the  element  changes  to  the  burning  state  by  the  propagation  of  the  fire, 
or  (3)  the  element  is  charred.  Once  the  material  begins  flaming  combustion, 
it  is  assumed  to  emit  smoke,  toxic  gases,  and  heat  until  the  material  is 
charred.  Thus  the  model  assumes  that  the  flaming  state  is  self-sustaining 
while  the  smoldering  state  requires  a  minimum  level  of  exterior  heating  to 
continue. 


2.2  MODELING  THE  CABIN  ATMOSPHERE 

DACFIR3  uses  the  "zone"  method  of  modeling  the  accumulation  of  com¬ 
bustion  products  in  the  cabin  atmospheres.  Zone  models  have  been  used  by  a 
number  of  investigators  for  analyzing  fires  in  residential  buildings  in  the 
past  few  years  and  this  method  is  developing  into  somewhat  of  a  standard  tech¬ 
nique  (reference  11).  The  method  has  proven  to  be  a  practical  first  approxi¬ 
mation  for  the  complex  problem  of  the  fluid  dynamics  and  heat  transfer  of  fire 
in  room-sized  enclosures. 

The  basic  assumption  of  zone  modeling  Is  that  the  hot  combustion  pro¬ 
ducts  rising  from  the  flames  gather  at  the  compartment  ceiling  to  form  a  well- 
mlxed  zone  or  layer.  As  the  fire  progresses,  more  and  more  of  the  original 
uncontaminated  gas  In  the  lower  region  of  the  compartment  Is  "pumped"  through 
the  flames  and  plume  to  enter  the  hot  upper  zone.  As  long  as  the  compartment 
volume  is  not  large,  the  turbulent  mixing  caused  by  the  rapidly  rising  fire 
gases  Is  sufficient  to  keep  the  upper  products  zone  well -mixed,  which  results 
In  a  uniform  temperature  and  composition.  This  model,  generally  supported  by 
experimental  observations,  allows  a  simple,  lumped  parameter  treatment  of  the 
conservation  equations  for  the  cabin  atmosphere.  Single  values  of  tem¬ 
perature,  density,  and  species  mass  fractions  are  used  to  describe  the  upper 
and  lower  zones.  The  boundary  between  the  zones,  the  "thermal  discontinuity," 
is  assumed  to  be  a  horizontal  plane  parallel  to  the  cabin  floor.  The  volumes 
of  the  upper  and  lower  zones  are  defined  by  the  position  of  this  plane  and  the 
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walls,  ceiling,  floor,  and  partitions  of  the  compartment.  Figure  2  shows  the 
main  features  of  the  two-zone  approximation.  When  the  cabin  Is  divided  into 
several  compartments,  the  two- zone  approximation  applies  to  each  compartment— 
when,  of  course,  conditions  are  such  that  some  hot  gases  have  penetrated  into 
every  compartment. 


- '  I  Upper  Zone 
Thermal  _ . 

Discontinuity  Entrainment 
Lower  Zone 


Flow 
to/ from 
Adjacent 
Compartment 


"77777777 


7777/777 


T  ,  T  *  upper  and  lower 
zone  temperatures 

xiu*  xit  =  upper  and  lower 
zone  concentra¬ 
tion  of  the  jth 
gas  specie 


Figure  2.  Zone  Model  of  the  Compartment  Atmosphere. 


The  simulation  of  the  behavior  of  the  cabin  atmosphere  consists  of 
numerically  solving  the  conservation  equations  for  each  zone  in  one  or  more 
compartments  for  the  values  of  temperature,  smoke  and  species  concentration, 
and  other  gas  properties  as  simulated  time  progresses.  Heat,  smoke,  and  gas 
accumulation  are  of  primary  interst  because  of  their  effect  on  emergency 
evacuation.  We  are  also  Interested  In  the  thermal  radiation  from  the  hot 
upper  zone  to  the  burning  or  about-to-be-lgnlted  materials.  This  radiation 
will  supplement  the  local  flame  radiation  to  Increase  the  flame  spread  and 
emission  rates  of  the  materials.  In  this  way,  the  cabin  atmosphere  Is  coupled 
to  the  fire  development  processes.  Details  of  the  cabin  atmosphere  model  are 
given  in  Section  4. 
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TABLE  1 

MATERIAL  FLAMMABILITY  CHARACTERISTICS  USED  AS  INPUT  FOR  DACFIR3 


Symbol 


Description 


fh  Horizontal  flame  spread  rate, 

f  Upward  flame  spread  rate, 

f .  Downward  flame  spread  rate. 

tf  Time  interval  from  flame  contact  to  the  start  of  flaming 

combustion. 

t^  Time  interval  from  the  start  to  the  end  of  flaming 

combustion. 

r^  Heat  release  rate  per  unit  area  during  flaming. 

r  f  Smoke  release  rate  per  unit  during  flaming. 

rf(i)  Release  rate  per  unit  area  of  the  ith  gas  specie  during 

flaming. 

qp  Heat  flux  at  which  the  material  begins  to  smolder. 

t  '  Time  interval  from  the  time  at  which  the  exosure  flux 

p  reaches  qp  to  the  start  of  smoldering. 

t  Time  Interval  from  the  start  of  smoldering  to  the 

p  consumption  of  the  material. 

t  Time  required  for  smoldering  to  stop  after  the  heat  flux 

p  falls  below  qp. 

r$s  Smoke  release  rate  per  unit  area  during  smoldering. 

r$(1)  Release  rate  of  the  1th  gas  specie  dw.’ing  smoldering. 

AHc  Effective  heat  of  combustion. 

y  Stoichiometric  oxygen-to-fuel  mass  ratio 

0^  Pyrolyzate  vapor  density. 

Pyrolyzate  vapor  surface  blowing  velocity. 

R.  Radiated  fraction  of  AH 

f  c 


2.3  RADIATION  MODELING 


As  explained  above,  the  primary  parameter  determining  the  fire  growth 
rate  in  DACFIR  is  the  level  of  thermal  radiation  reaching  the  burning  and 
smoldering  materials.  Calculation  of  the  radiation  intensities,  however,  is  a 
difficult  problem.  While  much  work  has  been  done  on  the  radiation  of  medium- 
scale  diffusion  flames,  the  results  are  still  limited  to  simple  geometries, 
e.g.,  horizontal  pool  fires  or  vertical  wall  fires.  Added  to  this  is  the  fact 
that  precise  radiation  calculations  are  inherently  complicated  and  highly 
dependent  on  the  exact  flame  geometry  and  temperature.  For  these  reasons  the 
DACFIR  model  must  use  simple,  approximate  relationships  for  flame  radiation. 
The  computations  assume  a  cylindrical  flame  volume  of  uniform  fixed 
temperature;  spectrally  gray  emissivity;  and  a  height  proportional  to  the 
hydraulic  radius  of  the  base  area.  Further,  for  materials  within  the  base 
reagion,  no  variation  of  flux  with  position  is  assumed;  and  an  average  value 
for  all  elements  in  the  base  is  used.  This  method  cannot  be  regarded  as  pre¬ 
cise  in  any  detail,  but  it  provides  an  indication  of  the  relative  radiation 
intensity  and  therefore  of  the  materials'  behavior.  This  level  of  approxima¬ 
tion  is  justified  partly  by  the  practical  problems  encountered  when  more 
detailed  methods  are  used  and  partly  by  the  degree  to  which  the  materials' 
properties  are  known  as  functions  of  radiation. 


2.4  ORGANIZATION  OF  THE  COMPUTATIONS 

The  DACFIR  model  is  Implemented  as  FORTRAN  computer  program  wnich 
integrates  the  differential  equations  and  simulates  the  fire  spread  for  a 
given  set  of  initial  conditions.  The  solution  of  the  differential  equation 
set  describing  the  cabin  atmopshere  requires  an  implicit  integration  method 
and  a  relatively  small  time  step,  on  the  order  of  one  second  or  less.  The 
simulation  of  fire  growth,  however,  proceeds  somewhat  more  slowly— for  most 
cases,  several  seconds  may  pass  between  the  Ignition  of  new  material  elements. 
The  program  is  designed  to  deal  efficiently  with  these  different  time  scales 
by  using  a  nested  loop  structure  as  shown  in  Figure  3.  During  a  run  integra¬ 
tion  of  the  cabin  atmosphere  equations  proceeds  at  a  given  time  step.  At.  At 
a  pre-set  Interval,  usually  five  seconds,  the  program  reviews  the  status  of 
all  burning  and  smoldering  elements  and  makes  the  required  ignitions, 
burn-outs,  changes  of  emission  rate,  etc.  This  procedure  can  continue  until 
all  the  cabin  materials  considered  by  the  model  have  burned-up,  but  normally 
we  are  only  Interested  in  about  the  first  ten  minutes  of  the  fire. 

The  other  Important  point  shown  In  Figure  3  is  the  loop  structure  for 
dealing  wtlh  Individual  fires.  Subroutine  SCAN  identifies  connected  groups  of 
burning  elements  which  form  the  base  planes  of  fires.  As  each  group  is  iden¬ 
tified  and  isolated,  the  flame  volume,  local  radiation  levels,  and  the  spread 
of  burning  caused  by  this  group  are  found  by  the  subroutines  which  follow 
SCAN.  If  there  are  additional  burning  elements,  control  returns  to  SCAN  to 
build  the  additional  fire  bases  and  compute  the  spread  caused  by  them.  When 
all  elements  on  all  surfaces  have  been  searched  and  the  new  growth  computed, 
the  new  total  rates  of  heat  and  mass  emission  are  updated  prior  to  entering 
the  cabin  atmosphere  routines  controlled  by  subroutine  ATMOS. 
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Figure  3.  Organization  of  the  DACFIR3  Program. 
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SECTION  3 

FIRE  SPREAD  AND  PRODUCT  EMISSION 


The  spread  of  burning  and  smoldering  over  the  cabin  materials  and  the 
rates  at  which  combustion  products  are  emitted  are  modeled  in  0ACFIR3  using  a 
discrete  states-transition  method.  This  section  describes  the  assumption  and 
oepration  of  this  method. 


3. 1  ELEMENT  STATES 

Each  material  element  in  the  cabin  Is  considered  to  be  in  one  of  seven 
states.  The  state  of  an  element  may  change  with  time  depending  upon  the  size 
and  location  of  the  fire  or  fires  in  the  cabin  and  upon  the  properties  of  the 
material.  The  seven  element  states  are: 


State  1  (Virgin) 

The  material  comprising  the  element  is  in  its  virgin  state;  that 
is,  it  has  not  been  directly  affected  by  the  fire  or  fires. 


State  2  (Smoldering) 


The  material  comprising  the  element  is  smoldering;  (undergoing 
nonflaming  decomposition). 

State  3  (Flaming) 

The  material  comprising  the  element  is  undergoing  flaming 
combustion. 

State  4  (Charred) 

The  material  comprising  the  element  has  burned  out  and  will  no 
longer  smolder  or  burn. 


State  5  (Heating,  no  flame  contact) 

The  material  comprising  the  element  is  receiving  heat  flux  suf¬ 
ficient  to  cause  It  to  smolder  but  because  of  the  element's  thermal 
inertia  smoldering  has  not  yet  begun. 


State  6  (Heating,  with  flame  contact) 


The  material  comprising  the  element  is  being  touched  by  the  flames 
of  a  fire  but  has  not  yet  Ignited  because  of  the  element's  thermal 
inertia. 


State  7  (Smoldering  and  cooling) 

The  material  element  began  smoldering  when  the  heat  flux  to  the 
element  reached  a  specified  level;  the  heat  flux  has  now  dropped 
below  that  level  but  the  material  is  still  smoldering.  The  element 
will  continue  to  smolder  for  a  specified  time  and  then  transform  to 
the  charred  state  (State  4). 

There  are  two  emission  rates  associated  with  each  material,  a  smoldering 
rate  and  a  flaming  rate.  The  elements  in  States  2  and  7  emit  smoke  and  toxic 
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gases  at  the  smoldering  rates.  The  elements  in  State  3  emit  smoke,  heat,  and 
toxic  gases  at  the  flaming  rate. 

Element  State  1  is  a  beginning  state  in  which  no  transitions  have  yet 
occurred.  State  4  is  a  final  state  from  which  no  additional  state  transitions 
can  occur.  States  5  and  6  are  intermediate  states  in  the  transitions  to 
States  2  and  3.  The  relationship  among  the  states  is  given  in  the  transition 
diagram  of  Figure  4.  The  criteria  governing  state  transitions  are  described 
in  the  next  section. 


3.2  ELEMENT  STATE  TRANSITIONS 

An  element's  transition  from  one  state  to  another  is  governed  by  the 
properties  of  the  material  and  by  the  element's  relationship  to  the  fire  or 
fires  in  the  cabin.  A  fire  is  defined  as  a  set  of  contiguous  elements  in 
State  3.  The  fire  is  assumed  to  have  a  flame  volume  whose  size  and  radiation 
are  a  function  of  the  area  of  the  fire  base  and  various  properties  of  the  fuel 
and  the  cabin  atmosphere.  At  any  time  there  are  four  sets  of  elements  that 
are  candidates  for  state  transition.  They  are: 

1.  Nonburning  elements  adjacent  to  elements  that  are  flaming.  These 
elements  are  candidates  for  immediate  ignition. 

2.  Elements  that  are  inside  the  flame  volume  of  a  fire.  All  such  ele¬ 
ments  are  candidates  for  immediate  ignition. 

3.  Elements  in  the  general  vicinity  of  a  fire.  These  elements  are  not 
in  danger  of  immediate  ignition  but  are  candidates  for  transition 
to  smoldering  due  to  the  high  local  radiation  flux. 

4.  Elements  that  are  smoldering  or  flaming.  These  elements  are  can¬ 
didates  for  transitions  to  the  charred  state. 

The  criteria  for  the  state  transitions  that  can  occur  for  the  elements  within 
each  of  thse  sets  are  described  below. 

3.2.1  Ignition  by  Creeping  Flame  Spread 

The  rate  at  which  a  flame  front  propagates  over  a  surface 
depends  upon  several  factors:  the  type  of  material  at  the  flame  edge,  the 
size  of  the  flames,  the  orientation  of  the  surface,  and  the  external  radiation 
level.  The  flame  spread  rates  for  a  given  material  are  supplied  to  the  model 
in  tabular  form  as  functions  of  radiation.  The  flux  to  elements  adjacent  to 
flaming  elements  is  the  sum  of  the  radiation  from  the  adjacent  fire  and  the 
radiation  from  the  hot  gases  in  the  upper  zone.  Three  flame  spread  rates  are 
associated  with  a  vertical  surface:  vertical  upward  spread  (fu),  vertical 
downward  (f^)  and  horizontal  (f ^ ) .  One  flame  spread  rate  (f^,)  is  associated 
with  horizontal  surfaces.  The  rates  and  directions  are  shown  in  Figure  5. 

The  state  of  an  element  adjacent  to  a  flaming  element  will  be 
changed  to  State  3  after  the  original  flaming  element  has  been  in  State  3  for 
a  time  interval  equal  to  or  greater  than  d/fj  where  d  is  the  distance  between 
the  center  of  one  element  and  the  center  of  an  adjacent  element  (six  inches) 
and  f<  is  the  flame  spead  rate  in  the  appropriate  direction  for  the  flaming 
material  for  a  given  heat  flux.  Thus  the  flame  front  progresses  in  steps  of 
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Figure  4.  Element  States  and  State  Transitions. 
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Figure  5.  Flame  Spread  Routes  for  Elements  on 
Vertical  and  Horizontal  Surfaces. 


six  inches  at  a  rate  determined  by  the  material  characteristics  and  the 
radiation  environment.  By  this  process--cal led  "creeping  flame  spread"  to 
distinguish  it  from  fire  jumping  across  nonconnected  surfaces— a  fire  can  pro¬ 
pagate  around  the  interior  lining  of  the  cabin  or  over  a  seat.  The  fire  may 
also  propagate  by  this  process  from  a  seat  adjacent  to  a  sidewall  to  the 
sidewall  and  vice  versa.  A  special  distance  and  rate  for  this  spread  are 
supplied  in  the  program  input.  If  a  given  material  will  not  burn,  the  flame 
spread  rate  can  be  given  as  zero  and  the  fire  will  not  propagate  over  the 
material.  The  fire  also  will  not  propagate  to  any  element  in  State  4 
(charred).  The  state  transitions  which  can  occur  by  creeping  spread  are  all 
transitions  to  State  3  from  any  of  States  1,  2,  5,  or  6.  These  state  tran¬ 
sitions  are  shown  in  Figure  4. 

The  flame  spread  rates  used  to  determine  if  a  state  transition 
occurs  are  functions  of  the  radiation  flux  arriving  at  the  material  at  the 
edge  of  the  fire.  This  flux  is  designated  as  q^  (Btu/ft^-sec)  and  is  computed 
as  a  function  of  the  flame  volume,  emissivity,  and  emissive  power,  and  of  the 
radiation  from  the  upper  gas  zone.  The  details  of  the  computation  of  q^  are 
given  in  Section  5. 

3.2.2  Ignition  by  Flame  Contact 

The. flames  of  a  fire  on  one  surface  may  touch  other  surfaces 
that  are  not  physically  connected  to  the  surface  containing  the  fire  base. 
Elements  on  the  surface  In  contact  with  the  flames  will  eventually  ignite, 
given  sufficient  exposure.  In  OACFIR,  fires  can  propagate  by  flame  contact  to 
any  surface  directly  above  a  burning  element,  provided  that  the  current  flame 
length  makes  contact  possible.  For  fires  originating  on  the  vertical 
sidewalls,  flames  may  contact  and  ignite  the  first  row  of  elements  on  the 
celling  or  underside  of  a  hatrack  or  passenger  service  unit  (PSU),  but  no 
extension  of  the  flames  horizontally  beneath  the  surface  is  considered. 

The  equation  for  computing  the  flame  length  is  presented  in 
Section  5.  For  a  fire  on  the  floor  the  height  is,  of  course,  measured  from 
the  floor.  For  fire  on  the  sidewall  the  height  is  measured  from  the  midpoint 
between  the  upper  and  lower  extremes  of  the  fire.  On  a  seat,  the  height  is 
measured  from  a  point  which  is  weighted  average  of  the  midpoints  of  four  sur¬ 
faces  which  make  up  the  seat  backrest.  The  weight  for  each  surface  is  the 
number  of  elements  on  that  surface  in  State  3  divided  by  the  total  number  of 
elements  in  State  3  on  the  backrest. 

The  state  transitions  which  occur  for  elements  in  contact  with 

flames  are: 

•  Elements  in  State  1  go  to  State  6 

•  Elements  in  State  2  go  to  State  3 

•  Elements  in  State  5  go  to  State  6 

•  Elements  In  State  6  go  to  State  3,  i_f  the  element  has  been 

In  State  6  longer  than  a  time  interval  tf. 

The  time  tf,  the  ignition  delay  for  flame  contact  ignition,  is  specified  in 
the  Input  for  each  material  as  tabular  function  of  heat  flux.  The  radiation 
level  qg,  experienced  by  elements  within  a  flame  volume  is  a  function  of  the 
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same  quantities  as  qj  but  will  be  of  greater  magnitude.  The  modeling  of  q£  is 
also  discussed  in  Section  5. 

3.2.3  Transition  to  Smoldering 

Elements  in  the  vicinity  of  a  fire  may  receive  sufficient  heat 
flux  from  the  flames  to  begin  smoldering.  In  DACFIR3,  this  process  is  allowed 
for  elements  on  the  floor  In  the  vicinity  of  a  fire  on  the  floor,  elements  on 
the  top  of  a  seat  cushion  in  the  vicinity  of  a  fire  on  the  seat  cushion,  ele¬ 
ments  on  the  bottom  of  a  seat  cushion  above  a  fire  on  the  floor,  and  elements 
on  the  ceiling  surfaces  above  a  fire  on  a  seat.  Associated  with  each  material 
is  an  input  parameter,  qp,  which  defines  the  heat  flux  level  above  which  the 
material  will  begin  to  smolder  after  a  brief  radiation  exposure. 

Elements  on  the  floor  In  the  vicinity  of  a  fire  on  the  floor  and 
elements  of  a  seat  cushion  In  the  vicinity  of  a  fire  on  the  cushion  top  will 
receive  a  heat  flux  level  greater  than  or  equal  to  qp  if  they  are  within  a 
distance  xp  of  the  fire  base  edge.  This  distance  is  referred  to  as  the 
"smoldering  range"  of  the  fire,  and  is  found  as  a  function  of  the  fire 
radiation.  The  radiation  flux  in  this  case  is  designated  as  q3.  Expressions 
for  q3  and  xp  are  given  in  Section  5. 

The  following  state  transitions  occur  for  elements  considered  to 
be  within  the  smoldering  range  of  the  fire: 

•  Elements  In  State  1  are  changed  to  State  5 

•  Elements  in  State  6  are  changed  to  State  5 

•  Elements  in  State  5  are  changed  to  State  2  j_f  they  have  been 

In  State  5  for  a  period  greater  than  tp. 

The  time  tp  Is  an  Input  parameter  and  is  a  function  of  the  type  of  material 
comprising  the  element. 

3.2.4  Transitions  of  Elements  in  the  Smoldering  or  Flaming  States 

For  each  time  step  that  an  element  is  in  State  2  (smoldering)  or 
State  3  (flaming)  a  certain  percentage  of  the  total  combustion  products  that 
can  be  emitted  are  emitted.  For  elements  in  State  2  this  percentage  is  calcu¬ 
lated  by  the  expression 

P2  »  (At/tpc)  100  (3-1) 

where  At  is  the  length  of  time  step  and  tpc  is  an  input  parameter  defining  the 
total  time  required  for  the  element  to  become  charred  at  a  heat  flux  of  qp. 

For  elements  In  State  3,  the  percentage  Is  calculated  by 

P3  *  (At/tfc)  100  (3-2) 

where  tfr  is  an  input  parameter  defining  the  total  time  required  for  the 
element  to  become  charred  due  to  flaming  combustion.  The  parameters  tpc  and 
tfc  are  functions  of  the  material  type  and  heat  flux.  The  heat  flux  received 
by  a  flaming  element  Is  q2- 

Each  time  P?  and  P3  are  calculated  for  an  element  they  are  added  to 
the  total  percentage  of  decomposition,  Pj,  for  the  element.  The  value  of 
Pt  Is  initially  zero  and  increases  by  P?  or  P3  for  each  time  period  that  the 
element  is  in  State  2  or  State  3.  When  Pj  for  an  element  becomes  equal  to  or 
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greater  than  100,  the  element  state  Is  changed  from  State  2  or  State  3  to 
State  4,  charred. 

While  an  element  Is  In  State  2,  It  emits  smoke  and  toxic  gases  at 
rates  which  are  Input  to  the  program  for  each  material.  While  an  element  Is 
In  State  3,  It  emits  heat,  smoke,  and  toxic  gases  at  rates  which  are  deter¬ 
mined  from  an  Input  table  of  these  rates  versus  heat  flux,  <\£. 

Elements  that  are  In  State  2  reach  this  state  because  the  heat  flux 
being  received  by  the  element  Is  greater  than  qp.  If  an  element  In  State  2  Is 
found  at  some  later  time  to  be  outside  the  smoldering  range  of  a  fire.  It  will 
no  longer  receive  a  heat  flux  greater  than  qp.  In  this  case  the  element  state 
Is  changed  from  State  2  to  State  7.  While  In  State  7  the  element  continues  to 
smolder  until  it  has  been  In  this  state  for  a  preset  time  tpe.  The  time 
tpe  is  an  Input  parameter  for  each  material.  Once  the  element  has  been  In 
State  7  for  a  time  tpe,  It  stops  smoldering  and  Its  stat*  Is  changed  from 
State  7  to  State  4.  *rhe  change  to  State  4  occurs  even  though  the  element  Is 
not  necessarily  completely  consumed.  Nonetheless,  no  ot  .er  transitions  are 
allowed  for  the  element.  Since  the  time  period  for  the  simulation  Is  normally 
not  long  enough  for  the  fire  to  recede  from  a  smoldering  element  and  lthen  to 
grow  back  to  It  again,  this  approximation  should  be  adequate  for  most  pur¬ 
poses. 

One  other  state  transition  Is  possible.  An  element  that  Is  within 
the  smoldering  range  of  a  fire  but  has  not  begun  to  smolder  (because  a  time 
Interval  tp  has  not  elapsed)  Is  In  State  5.  If,  due  to  some  change  In  the 
fire  condition,  an  element  In  State  5  Is  found  to  no  longer  be  within  the 
smoldering  range  of  a  fire.  Its  state  Is  reset  to  State  1. 


All  the  state  transitions  that  are  possible  and  the  characteristic 
time  for  each  transition  have  been  shown  in  Figure  4.  The  transitions  depend 
upon  the  specified  material  parameters:  three  flame  spread  rates  (f^,  fu, 
f{l),  five  characteristic  times  (tf,  tfr,  tp,  tpc,  tpe),  and  a  characteristic 
heat  flux  level  (qp).  The  flame  spread  rates  and  the  first  two  time  constants 
are  taken  to  be  functions  of  heat  flux  as  given  In  the  program  Input. 


3.3  DETAILS  OF  THE  CABIN  GEOMETRY 

The  cabin  of  a  commerlcal  transport  airplane  may  range  In  size  from 
about  10  feet  In  width  by  50  feet  In  length  to  over  20  feet  In  width  and  more 
than  20  feet  In  length.  The  number  of  seats  may  range  from  about  90  to  over 
500.  .  In  the  larger  wlde-body  transports,  the  cabin  Is  divided  into  smaller 
compartments  by  lavatories  and  galleys.  For  even  the  smallest  cabin,  the 
total  material  surface  area  (seats  and  cabin  lining  surfaces)  exceeds  5000 
square  feet.  All  the  surfaces  of  a  cabin  could  be  Included  In  the  grid  scheme 
of  0ACFIR3  but  most  of  the  elements  would  never  enter  Into  the  calculations. 
Therefore,  to  decrease  computer  time  and  storage  requirements,  DACFIR3  con¬ 
siders  only  a  three-seat  row-long  section  of  the  cabin  for  the  fire  spread 
simulation.  In  the  time  period  of  Interest  It  Is  not  likely  that  a  survlvable 
fire  will  spread  beyond  this  region. 

To  describe  the  cabin  in  the  Input,  the  user  of  DACFIR3  need  only 
supply  a  few  overall  dimensions  of  the  cabin  lining  surfaces  (i.e.,  the  floor. 
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sidewalls,  hatracks/stowbins,  etc.)  and  the  components  of  the  unit  vector  nor¬ 
mal  to  the  surface.  The  program  assembles  the  surfaces  in  the  proper 
arrangement,  divides  the  surfaces  Into  the  unit  elements,  and  sets  all  indices 
and  counters  necessary  for  the  fire  simulation  to  proceed  over  the  surfaces. 
All  common  transport  category  cabin  Interior  linings  can  be  represented  by 
this  method. 

All  dimensions,  locations,  and  directions  in  DACFIR3  are  specified  in  a 
single  coordinate  system.  A  right-handed  cartesian  system  is  used  with  the 
origin  located  in  the  forward  lower  right-hand  corner  of  the  cabin  as  viewed 
facing  forward.  This  coordinate  system  Is  shown  in  Figure  6.  The  figure  also 
shows  a  typical  location  for  the  "detailed  section,"  the  part  of  the  cabin 
where  the  grid  scheme  is  applied. 

3.3.1  Cabin  Lining  Surfaces 

The  lining  surfaces  in  an  aircraft  cabin  include  the  carpet, 
sidewalls,  window  reveals  and  transparencies,  passenger  service  units,  stowage 
bins,  and  ceiling  panels.  Older  aircraft  may  use  the  shelf-like  hatracks 
Instead  of  the  stowage  bins.  Lining  surfaces  are  assumed  to  be  oriented 
either  vertically  or  horizontally  with  respect  to  the  floor.  A  maximum  of  20 
individual  lining  surfaces  can  be  used.  The  vector  normal  to  a  surface  is 
specified  to  identify  the  surface  orientation  and  a  "z  displacement"  is  used 
to  locate  the  bottom  edge  of  the  surface  with  respect  to  the  cabin  floor.  For 
horizontal  surfaces,  the  z  displacement  Is  the  distance  from  the  floor  to  any 
point  on  the  surface  (as  all  points  are  equidistant  from  the  floor).  A  sur¬ 
face  width  is  specified  which  is  the  surface  z  dimension  for  vertical  surfaces 
and  the  x  dimension  for  horizontal  surfaces.  The  cabin  cross  section 
(parallel  to  the  x-z  plane)  is  assumed  to  be  the  same  for  all  compartments. 

3.3.2  Seats 

Seat  positions  and  seat  row  widths  are  user-definable  in 
DACFIR3.  Seats  are  modeled  by  a  simple  bench-like  approximation.  Seats  side 
by  side  In  a  row  are  regarded  as  a  single  bench  seat,  called  a  seat  group, 
whose  width  may  be  specified  provided  that  it  Is  an  integer  multiple  of  0.5 
feet.  The  other  dimensions  of  a  seat  are  assumed  to  be  constant.  The  seat 
cushion  Is  1.5  feet  from  back  to  front;  the  backrest  3.5  feet  high,  and  both 
cushion  and  backrest  are  0.5  feet  thick.  Figure  6  shows  the  seat  con¬ 
figuration  and  dimensions.  Up  to  nine  seat  groups  may  be  represented  in  the 
model,  each  with  a  different  width  If  required.  All  seats  are  assumed  to  face 
forward,  with  the  vertical  backrest  planes  perpendicular  to,  and  the  horizon¬ 
tal  cushion  planes  parallel  to  the  floor.  The  location  of  each  seat  group  Is 
given  by  specifying  the  x  and  y  coordinates  of  the  forward  right-hand  corner 
of  the  seat  cushion  as  viewed  by  a  person  sitting  in  the  seat. 
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SECTION  4 

CABIN  ATMOSPHERE  MODEL 


The  objective  of  the  cabin  atmosphere  model  of  DACFIR3  is  to  predict  the 
changes  in  gas  temperature  and  composition  as  the  fire  develops.  The 
"two-zone"  method  of  representing  the  cabin  atmosphere,  used  in  the  earlier 
versions  of  DACFIR,  is  extended  in  DACFIR3  to  cabins  divided  into  several 
compartments.  The  two-zone  method  assumes  the  cabin  atmosphere  within  a  com¬ 
partment  to  be  divided  into  two  separate  horizontal  regions  or  zones.  The 
upper  zone  consists  of  combustion  products  and  heated  air  which  have  risen 
toward  the  ceiling  by  natural  buoyancy.  Below  this  region  is  the  lower  zone 
consisting  of  cooler,  denser,  less  contaminated  air  that  was  originally  in  the 
compartment  or  which  enters  the  compartment  through  doors  or  vents.  Both 
zones  are  assumed  to  be  well  mixed  so  that  single  values  of  temperature, 
smoke,  and  gas  concentrations  suffice  to  describe  the  state  of  each  zone. 
Changes  in  the  mass,  composition,  energy,  and  volume  of  each  zone  occur  by  the 
action  of  fires  in  the  compartment,  by  radiative  and  convective  heat  exchanges 
between  the  zones  and  with  the  confining  walls,  and  by  the  effects  of  flows  in 
and  out  of  vents. 

This  section  presents  the  governing  equations  for  the  DACFIR3  model  of 
the  cabin  atmosphere.  Much  of  this  development  is  similar  to  that  previously 
reported  for  the  earlier  versions  of  the  model  and  by  others.  Enough  new 
material  has  been  added,  however,  to  justify  a  complete  restatement  of  the 
physics  and  mathematics. 


4.1  ZONE  MODEL  FOR  MULTIPLE  COMPARTMENTS 

The  extension  of  the  two-zone  model  to  multiple  compartments  is 
straightforward.  Each  compartment  is  divided  into  an  upper  and  lower  zone 
when  hot  combustion  gases  are  present.  Flows  between  compartments  are  driven 
by  density  and  pressure  differences  between  the  zones.  Figure  7  shows  the 
interaction  between  two  compartments  including  the  flow  from  one  of  the  com¬ 
partments  to  the  exterior.  The  broken  lines  in  the  figure  define  the  control 
volumes  for  each  zone.  When  the  burning  material  of  a  fire  base  is  located 
below  the  thermal  discontinuity  It  Is  convenient  to  extend  the  upper  zone 
control  volume  down  through  the  flames  and  plume  to  the  burning  surface.  Mass 
flow  across  the  boundaries  of  each  control  volume  occurs  at  the  planes  of 
vents  and  at  the  entraining  surfaces  of  flames  and  plumes. 
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Figure  7.  Gas  Zone  Control  Volumes  and  Mass  Flows 


We  make  the  following  simplifying  assumptions  about  the  compartmentation 
and  cabin  design: 

•  A  maximum  of  four  compartments  Is  allowed 

•  Compartments  are  arranged  linearly  along  the  fuselage  so  that  no  com¬ 
partment  is  connected  to  more  than  two  others 

•  The  floor  levels  in  all  compartments  are  equal 

•  The  cross-section  (in  the  x-z  plane)  of  all  compartments  is  identical 

•  Any  or  all  compartments  may  have  vents  (doors,  escape  hatches,  etc.) 
leading  to  the  outside  or  to  other  compartments,  but  no  compartment 
has  more  than  six  vents 

•  The  partitions  dividing  the  compartments  are  vertical  planes  of 
negligible  thickness. 

All  the  above  assumptions,  with  perhaps  the  exception  of  the  last,  are 
appropriate  to  any  modern  cabin  design.  We  also  assume  that  the  cabin  floor 
is  horizontal  so  that  the  plane  of  the  boundary  between  the  zones  is  always 
parallel  to  It. 

When  DACFIR3  Is  used  with  multiple  compartments,  the  user  specifies  the 
compartment  in  which  the  grid  scheme  for  tracking  interior  fires  is  used.  All 
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burning  and  smoldering  of  interior  materials  will  be  confined  to  this  compart¬ 
ment  and  the  limitations  of  the  grid.  The  activity  in  other  compartments  is 
limited  to  the  accumulation  of  combustion  products  by  flow  through  vents-- 
including  the  inflow  of  gases  from  an  exterior  fire--and  the  exchange  of  heat 
between  zones  and  to  the  confining  walls. 

The  equations  governing  the  behavior  of  each  zone  in  each  compartment 
are  presented  in  the  next  section. 

4.2  CONSERVATION  EQUATIONS 

For  the  control  volumes  defining  the  upper  and  lower  zones  in  a  compart¬ 
ment,  application  of  the  conservation  of  mass,  energy,  and  species  allows  us 
to  write  a  set  of  ordinary  differential  equations  for  the  mass,  composition, 
and  energy  in  each  zone.  When  we  add  several  additional  relationships  (the 
gas  law,  etc.)  and  apply  the  fact  that  the  volume  of  each  compartment  is 
fixed,  the  equation  set  is  closed  and  may  be  solved  numerically  for  specified 
initial  conditions. 

The  develoment  of  the  conservation  equations  follows  standard  procedure 
and  so  will  be  given  in  an  abbreviated  form. 

4.2.1  Conservation  of  Mass 

Conservation  of  mass  applied  to  both  zones  in  a  compartment 
yields,  for  the  upper  zone 


A.J  - 

dt  u 


l  Giu  +  l  Gn  +  l  Gf 

vents  plumes  K  fires 
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and  for  the  lower  zones 
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vents  plumes  K 
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The  superscript  i  indicates  the  ith  compartment  if  more  than  one  is  used.  The 
summations  of  the  various  mass  flow  rates  are  understood  to  be  taken  over  the 
surface  areas  of  each  control  volume  through  which  mass  passes.  Thus  the  ven< 

flows,  G^,  occur  only  if  some  of  the  zone  surface  is  in  contact  with  a  vent 

opening.  The  second  summation  on  the  right  of  both  equations  represents  the 
rate  of  entrainment  of  mass  into  the  plumes  of  fires  whose  bases  extend  into 
the  lower  zone.  The  assumed  algebraic  sign  for  all  mass  flows  and  thus  for 
Gp  is  positive  for  flow  into  the  control  volume.  Therefore  the  sum  is  posi¬ 
tive  in  (4-1)  and  negative  In  (4-2).  Vaporizing  fuel  adds  mass  to  the  upper 
zone  across  the  base  planes  of  any  burning  surface.  This  flow,  small  but  not 
always  negligible,  adds  to  the  right  hand  side  of  (4-1).  Further  analysis  of 
each  of  the  mass  flow  rates  will  be  given  in  later  sections. 

4.2.2  Conservation  of  Species 

In  DACFIR3  we  make  the  following  assumptions  about  the  com¬ 
position  of  the  cabin  atmosphere: 
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•  All  gas  consists  primarily  of  five  major  species  N2,  O2, 
"fuel  vapor,"  CO2,  and  H2O. 

•  Any  additional  species,  e.g.,  CO,  HCN,  HC1 ,  etc.  exist  at 
such  a  low  concentration  that  the  contribution  to  the  mean 
molecular  weight  in  any  zone  is  negligible.  These  are  the 
"trace  species." 

•  Smoke  is  treated  as  a  trace  specie  without  mass  but  with  a 
concentration  related  to  its  optical  density. 

These  assumptions  allow  a  unified  and  economic  treatment  for  the  computation 
of  smoke  and  gas  concentrations.  The  major  species  were  chosen  because  most 
cabin  materials  are  polymers  of  carbon,  hydrogen,  oxygen,  and  nitrogen  and  so 
might  be  expected  to  produce  primarily  CO2  and  H2O  upon  combustion.  The  spe¬ 
cific  trace  species  and  their  molecular  weights  are  defined  by  the  user  in 
connection  with  the  emission  rate  data  supplied  for  the  materials. 
Concentrations  of  the  trace  species  and  the  major  species  are  given  in  terms 

of  mass  fraction,  x1 »  of  the  jth  specie  in  each  zone: 
j 
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and  for  the  lower  zone. 
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The  subscript  j  refers  to  the  jth  specie. 

For  the  vent  flow  terms,  some  flows  may  be  out  of  a  zone  and  so 
will  use  the  mass  fraction  of  the  jth  specie  in  that  particular  zone.  In 
other  cases  the  flow  will  be  into  the  zone  from  an  adjoining  compartment  or 
from  the  exterior  so  that  the  mass  fraction  will  be  characteristic  of  the  zone 
where  the  flow  originates.  Zone  subscripts  and  compartment  superscripts  have 
therefore  not  been  placed  on  the  mass  fraction  variable  for  the  vent  flow 
terms. 


Flames  and  plumes  of  fires  in  the  lower  zone  entrain  the  species 
from  that  zone  and  so  Xjp  appears  on  the  right  of  both  equations.  The  final 

two  terms  of  (4-3)  are  the  mass  generation  rates  of  the  jth  species  assumed  to 
occur  at  the  base  planes  of  fires  and  over  the  areas  of  smoldering  materials. 
This  is  the  point  of  connection  between  the  product  generation  data  of  trace 
species  specified  for  the  cabin  materials  and  the  atmosphere  model.  The  sum¬ 
mations  for  the  generation  terms  are  developed  in  the  flame  spread  sections  of 
the  computer  code  and  delivered  to  the  cabin  atmosphere  as  a  single  term  which 
applies  only  to  the  compartment  containing  the  element  grid.  Although  the 
model  assumes  that  no  plumes  exist  above  smoldering  regions  in  the  lower  zone 
and  that  no  significant  mass  is  generated  by  smoldering',  we  place  all  the 
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trace  species  generated  by  smoldering  in  the  upper  zone,  since  it  is  likely 
that  they  will  rise  to  that  area  by  buoyancy. 

For  all  true  chemical  species,  the  variables  Xi  are  mass  frac- 

tions  and  thus  unitless.  For  smoke,  the  input  generation  data  is  in  units  of 
"particles  of  smoke"  generated  per  square  foot  per  second.  The  unit  "particle 
of  smoke,"  introduced  by  Smith  (Reference  12)  to  quantify  rates  of  smoke 
release,  is  defined  as  that  amount  of  smoke  which  when  confined  to  a  volume  of 
one  cubic  foot  reduces  the  transmission  over  a  one  foot  light  path  by  10  per¬ 
cent.  Units  for  xj  for  smoke  are  then  particles  per  lbm  of  mixture. 

In  the  case  of  oxygen  consumption,  the  term  Wq^  is  negative  and 

its  value  is  obtained  from  the  specified  values  of  the  heat  of  combustion  for 
each  material,  aHcx,  the  total  rate  of  heat  release  <3|<  over  all  the  burning 
area  of  this  material,  and  the  stochiometric  oxygen-to-fuel  mass  ratio,  Yfc, 

The  subscript  k  denotes  the  kth  material.  We  assume  an  overall  stoichiometric 
combustion  so  that 
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k  materials 


\  VAHck 
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where  the  units  of  Wg  are  pounds  of  Og  consumed  per  square  foot  per  second. 

Generation  rates  of  the  major  product  species  C02  and  H20  are 
computed  from  the  rate  of  mass  consumption  as  estimated  by  the  ratio  Qx/aHcx. 
With  a  specified  molecular  weight  for  each  material,  u>x,  the  rates  of  produc¬ 
tion  of  carbon  dioxide  and  water  vapor  are 


“io2  ■  i  <44'-\><iv*v  -  (,-6> 
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The  stoichiometry  assumed  in  this  process  is  that  mole  of  fuel,  of  whatever 
type,  yields  one  mole  each  of  C02  and  H20  upon  combustion.  [This  situation 
would  strictly  apply  only  If  all  materials  could  be  characterized  by  a  struc¬ 
ture  (CH2)_  so  that  the  overall  reaction  would  be  (CH2)n  +  (3n/2)02  +  nC02  + 
nH20.J 

In  DACFIR3  no  unburned  fuel  vapor  is  assumed  to  be  present  in 
either  zone  so  that  the  value  of  xfuel*  for  each  20ne  a^Wflys  zero. 
Similarly,  no  consumption  rate  term,  Wfue],  analogous  to  Wq^,  is  computed  and 

the  conservation  equation  for  the  fuel  species  Involves  only  null  arithmetic. 
This  approach  of  keeping  a  variable  and  equation  for  fuel  vapor  was  taken  in 
anticipation  of  later  refinement  to  include  such  phenomena  as  Incomplete 
combustion  and  the  accumulation  of  unburned  hydrocarbons  which  may  become 
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involved  in  other  fires.  In  the  same  sense,  there  are  no  production  or 
consumption  terms  for  the  lower  zone  species  equations  in  this  version. 
Consequently  mass  fractions  in  the  lower  zones  of  all  compartments  will  not 
change.  (This  assumes  all  lower  zones  exchange  flows  only  with  one  another. 
This  assumption  is  discussed  in  the  section  on  vent  flows.)  By  writing 
species  conservation  equations  for  the  lower  zone  we  anticipate  refinements  to 
Include  mixing  between  upper  and  lower  zones,  which  is  known  to  occur  but  for 
which  no  model  Is  yet  available  (see  Reference  13). 

4.2.3  Conservation  of  Energy 

In  applying  the  first  law  of  thermodynamics  to  each  zone  control 
volume  we  may  assume  that  the  heat  addition  or  loss  occurs  at  constant 
pressure.  There  are  no  work  Interactions  so  the  conservation  of  energy  for 
the  upper  zone  is 
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and  for  the  lower  zone 
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The  first  assumption  to  note  is  that  we  take  the  specific  heat  at  constant 
pressure,  cp,  to  be  Independent  of  temperature  and  composition.  We  use  the 
value  for  room  temperature  air,  0.24  Btu/lbm-R,  throughout.  As  discussed  for 
the  species  conservation  equations,  the  energy  transported  through  the  vents 
will  depend  on  the  flow  direction  so  that  the  gas  temperature  involved  must  be 
separately  determined  In  each  case.  When  the  equations  apply  to  the  compart¬ 
ment  containing  Interior  fire,  the  plume  entrainment  (if  any)  will  carry 
enthalpy  from  the  lower  zone  to  the  upper  as  given  by  the  second  term  on  the 
right  in  each  case.  The  third  term  on  the  right  of  (4-8)  represents  the  total 
heat  release  by  all  interior  fires  since,  by  convention,  all  fire  bases  contact 
the  upper  zone.  The  convective  loss  to  the  non-burning  surfaces  appears  as  an 
addition,  although  the  term  Is  computed  In  a  way  that  preserves  the  sense  of 
the  heat  flow  which  Is  normally  out  of  the  zone.  The  remaining  two  terms  are 

the  rates  of  radiation  absorbed  by  the  zone  gas,  $rin>  and  the  rates  of 

emission.  Grout*  The  method  of  computing  all  the  energy  flow  terms  is 
presented  In  the  following  sections. 
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4.3  INTERIOR  FIRES 

Fires  in  the  compartment  entrain  mass,  species,  and  energy  from  the 
lower  zone  into  the  upper  zone.  A  number  of  analyses  of  the  behavior  of  fires 
of  the  scale  considered  by  DACFIR  (i.e.,  base  diameters  on  the  order  of 
several  feet)  are  available  in  the  literature.  The  model  of  Steward 
(Reference  14)  and  its  its  extension  by  Fang  (Reference  15)  have  been  chosen 
for  use  in  DACFIR  because  it  incorporates  heat  generation  in  the  flames  and 
includes  the  effect  of  reduced  oxygen  concentration  in  the  entrained  gas. 

Steward's  model  was  developed  for  fires  with  horizontal,  circular  bases 
and  axi symmetric  flame  and  plumes  rising  above  the  base  plane.  DACFIR  allows 
fires  on  vertical  and  horizontal  downward  facing  surfaces  with  base  plane 
areas  of  any  pattern  that  can  be  formed  by  connected  square  elements.  Since 
we  employ  the  fire  model  only  to  estimate  the  entrainment  and  not  such  items 
as  flame  temperature  or  velocity  we  adopt  Steward's  entrainment  relationships 
with  the  understanding  that  better  models  which  include  all  the  complex 
geometric  effects  of  the  base  shape  and  orientation  are  needed.  For  the  fire 
base  radius,  y0,  in  the  Steward  model  we  use  the  hydraulic  radius  of  the  base 
defined  as 


yQ  =  2A/P  (4-10) 

where  A  is  the  base  area  and  P  is  the  base  perimeter. 

The  fire  structure  assumed  by  the  Steward  model  is  shown  in  Figure  8. 

An  axi symmetric  column  of  gas  rises  from  the  burning  surface,  the  column 
radius  varying  with  height  as  air  is  entrained  into  the  turbulent  flow.  Two 
regions  in  the  plume  are  Identified:  the  combustion  zone  which  extends  from 
the  base  to  the  point  at  which  a  stoichiometric  mixutre  is  reached,  and  above 
that,  the  plume  zone  where  no  further  reaction  is  assumed.  We  have  used  the 
term  plume  earlier  to  signify  the  entire  gas  column  that  is  a  fire,  but  here 
we  reserve  the  term  for  the  region  of  the  column  where  there  is  no  further 
heat  generation.  Note  here  that  in  Steward's  model  the  stoichiometric  mixture 
height,  Zs,  and  the  visible  flame  height,  hf,  do  not  coincide.  Further 
discussion  of  this  point  is  given  below. 


Figure  8.  Fire  Structure  and  Entrainment  Flows  in  Steward's  Fire  Model. 
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The  rate  of  mass  flow  into  the  combustion  zone  from  the  base  to  a  height 
z  is  given  in  Steward's  theory  by 

Gec  *  Auop^{[0.673  Ec4/5o>“3/5 ( 1-u.) 1/5 ( gy Q/u ^ ) 17 52/yQ  +  1]5/2  -  1}.  (4-11) 

where 


bl  = 


cvl/(vl +  xiofc 


)]  . 


The  values  of  u0,  y.  and  AHC  are  those  for  the  material  of  the  fire  base. 
The  combustion  zone  extends  from  the  fire  base  to  the  stoichiometric  mixing 
height  given  by 

Z$  =  1.49  Ec-4/5[u,/(l-a.)]1/5(<dPjlP0  +  T/xA02)2/5(Pouo/Pa/9yo)2/5  (4-12) 


Above  this  point  the  plume  entrains  lower  zone  gas  at  a  rate  proportional  to 
the  distance  above  zs  given  by 

Gep  =  W^Cl.09  E4£5  (gys/u2)1/5  (1-ps/p4)1/5  (z-zs)/ys  +  1]5/3  -  1} 

(4-13) 

where  ys,  us,  and  ps  are  the  column  radius,  velocity,  and  density  at  the 

stoichiometric  point.  As  is  the  column  cross-sectional  area  *y|.  The  values 
of  these  quantities  are  computed  from  expressions  given  by  Fang,  using  the 
same  information  about  the  fire  base  and  fuel  vapor  employed  in  (4-11,  12,  and 
13).  For  brevity  these  expressions  have  been  placed  in  Appendix  B  of 
Volume  II. 

DACFIR3  tests  for  the  postion  of  the  fire  base  midpoint  and  the  point  of 
stoichiometric  mixing  with  respect  to  the  thermal  discontinuity.  The 
appropriate  values  of  z  are  then  Inserted  in  (4-11)  and  (4-12  and  13)  as 
required  to  find  the  total  entrainment  by  this  fire  column.  The  same  calcula¬ 
tion  is  performed  for  each  interior  fire  whose  base  midpoint  lies  below  the 
thermal  discontinuity.  The  sum  over  all  the  fires  Is  the  mass  entrainment 
term  of  (4-1).  In  the  energy  and  species  equations,  the  mass  entrainment  rate 
is  multiplied  by  CpT\  or  the  species  mass  fractions  xlj  as  appropriate  to  find 
the  entrainment  flows  of  these  quantities. 

Two  added  points  should  be  mentioned  concerning  the  entrainment 
expressions.  First,  the  values  of  the  heats  of  combustion  used  in  (4-11) 
and  (4-12)  are  "effective"  values  wherein  the  standard  heat  of  combustion  is 
reduced  by  a  given  factor  to  account  for  the  heat  lost  by  flame  radiation. 

The  size  of  this  loss  is  not  well  defined  but  is  on  the  order  of  30  to  40  per¬ 
cent  for  the  smokey  flames  of  polymeric  materials.  Since  the  user  supplies 
the  effective  aHc  in  the  program  input,  no  specific  assumption  about  the 
radiation  adjustment  Is  made  by  DACFIR3.  Second,  the  values  of  the  entrain¬ 
ment  constants  Ec  and  Ep  are  not  well  defined.  Estimates  for  the  values  range 
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from  about  0.1  to  as  high  as  0.3.  The  values  chosen  for  DACFIR3  are  0.15  for 
Ec  and  0.10  for  Ep. 

The  Steward  model  is  also  used  to  compute  flame  length,  hf,  for  the  fire 
spread  computations.  The  flame  length  is  not,  as  one  might  first  expect,  the 
same  as  the  stoichiometric  height  z«  but  is  typically  several  times  this 
distance.  The  prediction  of  flame  length  was  developed  by  Steward  by  corre¬ 
lating  the  model  with  experiemental  data  or  fires  of  a  large  range,  of  sizes. 
Further  discussion  of  the  flame  length  prediction  method  is  given  in 
Section  5. 

4.4  VENT  FLOWS 

DACFIR3  has  two  methods  for  dealing  with  the  flow  of  gases  through  com¬ 
partment  vents.  In  the  first  method  the  user  specifies  the  volume  flow  rate 
and  direction  at  a  vent.  The  program  uses  this  constant  "forced"  flow 
throughout  the  simulation  regardless  of  the  compartment  atmosphere  conditions. 
Forced  flow  is  used  when  the  model  is  to  be  compared  to  the  results  of  tests 
which  employ  forced  ventilation.  It  is  also  used  in  the  special  case  of  an 
exterior  fire  at  the  vent.  In  this  case  the  temperature  and  composition  of 
the  inflowing  fire  gas  must  also  be  given  in  the  input. 

In  the  second  method,  the  program  computes  the  flow  rate  using  the  com¬ 
partment  pressure  and  the  buoyancy  of  the  hot  gas.  This  "free"  flow  situation 
is  that  which  would  normally  occur  at  cabin  exit  doors  and  between  compart¬ 
ments  in  a  post-crash  fire.  The  method  for  buoyancy  driven  flow  from  a  single 
compartment  to  the  exterior  single  compartment  to  the  exterior  is  now  a  stan¬ 
dard  technique  verified  by  a  number  of  experiments.  Extensions  of  the  method 
to  the  flow  between  two  compartments,  both  containing  layers  of  hot  gas,  were 
made  by  Tanaka  (Reference  16)  and  by  Emmons  (Reference  17).  The  extensions 
involve  no  new  physical  assumptions  but  require  a  complicated  switching  struc¬ 
ture  to  account  for  all  the  flow  situations  that  can  occur.  The  application 
of  the  free-flow  computations  to  multiple  compartments  in  DACFIR3  is  essen¬ 
tially  Emmons'  formulation  but  with  a  different  logical  structure  to  account 
for  all  possible  flow  cases. 

Figure  9(a)  shows  a  single  compartment  with  a  buoyant  flow  of  hot  gases 
to  the  outside  and  a  cold  return  flow  through  the  lower  part  of  the  vent.  In 
part  (b)  of  the  figure  the  interior  and  exterior  pressures  are  shown  as  func¬ 
tions  of  height  above  the  floor.  The  curves  are  given  by  the  hydrostatic  law 
for  constant  density  gas  in  each  zone.  The  rate  of  pressure  change  with 
height  in  the  interior  decreases  less  rapidly  above  the  thermal  discontinuity 
height,  z^,  since  the  density  in  the  upper  zone  is  less  than  that  in  the 
lower  zone.  In  the  example  of  Figure  9  (a)  and  (b),  the  lower  zone  tem¬ 
perature  and  density  are  equal  to  the  outside  values  so  that  the  slope  of  the 
pressure  curve  in  the  lower  zone  is  the  same  as  the  slope  outside.  The 
pressure  in  the  lower  zone,  however,  is  lower  at  every  height  than  the  outside 
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Flow 


(a)  &  (b)  Single  Compartment  Flow  to  Exterior 
(c)  &  (d)  Flow  Between  Two  Compartments 


Figure  9.  Pressure  and  Buoyancy  Driven  Flows  Through  Compartment  Vents. 
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pressure  because  of  the  entrainment  of  lower  zone  gas  by  fire  plumes  inside. 
This  results  in  a  lower  Interior  pressure  at  the  floor,  Pf,  than  the  outside 
pressure  at  the  same  level.  The  point  at  which  the  pressure  profiles  cross  Is 
called  the  neutral  axis,  zn,  which  defines  the  boundary  between  the  flows  in 
the  vent  plane. 

To  compute  the  flow  rates  for  this  situation,  the  Bernoulli  equation  is 
used  to  express  the  gas  velocity  at  the  vent  plane  In  terms  of  the  pressure 
difference  between  the  inside  and  outside.  The  product  of  the  velocity  and 
density  of  the  flowing  gas  is  the  mass  flow  rate  per  unit  area  as  a  function 
of  height.  Since  the  pressure  difference  is  a  piecewise  linear  function  of  z, 
analytic  forms  can  be  found  for  the  Integrals  of  the  mass  flow  rates  for  each 
section  of  the  pressure  difference  function.  The  flow  rates  are  corrected  in 
the  usual  way  by  use  of  an  empirical  orifice  coefficient.  The  mathematical 
details  of  this  procedure  are  given  below  in  the  specific  case  of  the  flow 
between  two  compartments. 

The  flow  between  two  compartments  can  result  in  a  more  complicated 
situation  as  shown  in  Figures  9(c)  and  (d).  Both  compartments  contain  two 
zones  and  the  pressure  ,t  the  floor  may  differ  between  compartments.  The  par¬ 
ticular  values  of  the  four  densities  involved  and  the  location  of  the  two-zone 
boundaries,  z^  and  z^?.  may  result  in  three  separate  neutral  axes,  zn»,  znK, 
and  znr  as  shown.  As  in  the  more  simple  case  of  9(a)  and  (b),  the  method  of 
computing  the  flow  rates  Involves  writing  the  mass  flux  as  a  function  of  the 
piecewise  linear  pressure-height  relationship.  Integration  over  z  from  the 
bottom  to  the  top  of  the  vent  gives  the  mass  flows.  In  the  case  shown  in  9(c) 
four  separate  flows  occur. 


The  mathematical  procedure  for  computing  the  flows  between  the  compart¬ 
ments  is  as  follows.  The  hydrostatic  law  applied  to  each  compartment 
(numbered  1  and  2)  gives  pressure  as  a  function  of  z  as 


Pl<2>  *  Pfl  - 

9 

2  i  zdi- 

(4- 14a) 
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Depending  on  the  values  of 

Pfl*  Pf2»  the  four 

densities,  and  the 

two  thermal 

discontinuity  heights,  one  to  three  neutral  axes  may  exist  or  there  may  be  no 
neutral  axes.  In  the  latter  case,  the  flow  will  either  be  all  into  or  all  out 
of  the  vent. 


To  find  the  neutral  axes,  if  any,  DACFIR3  solves  for  the  intersections 
of  each  of  the  lines  represented  by  equations  (4-14)  with  each  of  the  lines 
given  by  (4-15).  This  results  In  four  "candidate"  axis  heights  Zj ,  one  or 
more  of  which  Is  not  physically  possible.  The  case  where  the  upper-  or  lower- 
zone  densities  are  equal  In  both  compartments  and  thus  the  pressure  profiles 
are  parallel  Is  recognized  in  a  "dummy"  Intersection  height  of  a  negative 
value  Is  set.  The  set  of  four  candidate  axis  positions  is  tested  for  physi¬ 
cally  Incorrect  values— that  Is,  any  negative  value  (axis  below  floor  level) 
or  any  Intersection  on  the  wrong  side  of  the  thermal  discontinuity.  For 
example,  the  profile  of  the  upper  zone  pressure  in  compartment  1  could 
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intersect  the  profile  of  the  lower  zone  pressure  of  compartment  2  at  a  value 
of  z  less  than  z^,  but  since  the  upper  zone  of  1  does  exist  below  z^j  there 
is  not  a  neutral  axis  at  this  point. 

As  the  neutral  axes  are  identified  they  are  combined  in  a  list  with  the 
one  or  two  zVs.  When  the  list  is  complete  it  is  sorted  into  increasing 
numerical  order.  The  upper  and  lower  edges  of  the  vent  opening  are  then 
placed  in  the  list  and  any  or  zn  not  in  the  vent  opening  is  removed.  What 
remains  then  is  an  ordered  list  of  the  breaks  in  the  piecevise  linear  curve  of 
pressure  difference  between  the  compartments  versus  z  for  values  of  z  in  the 
vent  opening.  Symbolically,  this  function  is 

AP  =  fj(z) 

A P  -  f2(z) 

(4-16) 

• 

AP  *-  fn(z) 

The  functions  fj  are  either  constants  when  the  densities  on  either  side  of 
vent  are  equal  or  linear  functions  of  z  when  the  densities  are  not. 

For  the  case  where  AP  is  constant,  application  of  the  Bernoulli  equation 

gives 

Gk  *  /2  CAx(pAP)1/2  (zR  -  zk_1)  (4-17) 

where  G^  is  the  mass  flow  rate  through  the  vent  between  the  heights  z^i  and 
Zfc,  Ax  is  the  vent  width,  and  C  is  an  orifice  coefficient,  taken  as  0.6&  in  all 
cases.  The  value  of  p,  the  density  of  the  flowing  gas,  is  determined  from  the 
sign  of  AP  and  the  zone  from  which  the  gas  flows.  Equation  (4-17)  is  also  used 
in  the  special  case  when  AP  is  equal  to  zero. 

When  AP  is  a  function  of  z  the  mass  flow  rate  between  z^_i  and  z^  is 
Gk  *  (2/2/3)CAx  p1/Z  (l/a)[(APZk)3/2  -  (AP(zk_1))3/Z],  (4-18) 

where 


zo  <  2  <  zl* 

zl  <  2  <  z2* 


a  s  [AP(zk)  -  AP(zk.1)]/(zk  -  zk_1). 

As  In  (4-17)  the  specific  value  for  p  and  the  algebraic  sign  of  G^  depends 
upon  the  sign  of  the  pressure  difference  between  the  compartments.  The  deve¬ 
lopment  of  the  list  of  break  points  zu  ensures  that  AP(Z|()  and  aP(z|c_i)  have 
the  same  sign  (In  some  cases  one  of  the  values  may  be  zero)  so  that  the  choice 
of  flow  direction  and  density  Is  unambiguous. 

When  all  mass  flows  between  the  break  points  z^  have  been  found  they  are 
summed  Into  four  categories  depending  upon  the  origin  zone  and  destination 
zone  of  each  flow.  We  make  the  following  assumptions  about  the  behavior  of 
the  flows:  hot  gases  from  upper  zones  flow  only  to  adjacent  upper  zones  and 
cold  gases  from  lower  zones  flow  only  to  adjacent  lower  zones.  Between 
compartments  1  and  2  there  are  then  four  possible  net  flow  rates;  from  the 
upper  zone  of  1  to  that  of  2,  from  the  upper  zone  of  2  to  that  of  1,  and  the 
analogous  pair  for  the  lower  zones.  The  various  G^'s  are  summed  into  these 


30 


four  categories  and  the  corresponding  flows  of  energy  (heat)  and  species  are 
found  by  multiplying  the  G^'s  by  the  appropriate  values  of  cpT  and  xi  and 
summing.  p  J 

The  above  process  is  carried  out  for  every  vent  in  the  compartment  or 
compartments  at  each  time  step.  The  program  logic  is  constructed  so  that  no 
vent  need  be  considered  more  than  once. 

Two  special  situations  should  be  noted.  First,  when  the  vent  is  open  to 
the  outside  the  procedure  outlined  above  is  applied  by  assuming  the  exterior 
to  be  another  compartment  who^e  lower  zone  always  covers  the  vent.  The 
pressure  at  floor  level  in  this  "outside"  compartment  is  always  fixed  at  the 
user  specified  ambient  value.  Similarly,  the  outside  temperature  and  density 
are  set  to  fixed  values  in  the  input.  The  program  keeps  a  fixed  thermal 
discontinuity  height  for  the  outside  as  a  large  value,  much  greater  than  the 
top  of  any  vent. 

The  other  special  situation  is  the  case  of  an  exterior  fire  at  a  vent. 
Here  the  user  specifies  the  volume  Inflow  rate  of  exterior  fire  gas,  Vgxf; 
and  the  temperature,  density,  and  composition  of  the  gas,  T  Pexf,  and 
Xjexf*  These  values  are  fixed  throughout  the  run  irrespective  of  what 
conditions  develop  in  the  interior.  The  mass  inflow  rate  is  then  Ggxf  = 

pexf^exf  and  the  enep9y  and  sPecies  f1ws  are  cpTexfGexf  and  XjexfGexf* 

4.5  CONVECTION  AND  RADIATION 

Gases  in  both  the  upper  and  lower  zones  of  each  compartment  may  exchange 
heat  with  the  confining  surfaces  of  the  cabin  and  with  the  neighboring  zone. 
DACFIR3  incorporates  the  following  heat  transfer  mechanisms: 

•  Convective  heat  transfer  between  the  cabin  lining  surfaces  and  the 
gas  in  both  zones. 

•  Radiation  from  the  hot  upper  zone  gas  to  all  cabin  lining  materials. 

•  Radiation  to  each  zone  from  the  surfaces  in  contact  with  that  zone. 

•  Radiation  between  zones. 

The  other  possible  radiative  patterns  such  as  the  radiation  from  one  interior 
surface  to  another  are  assumed  to  be  Insignificant  and  are  ignored.  DACFIR3 
does  not  use  the  seat  surfaces  In  any  radiation  or  convective  heat  exchange. 

4.5.1  Convection 

The  convective  heat  transfer  Is  modeled  by  assuming  a  constant 
convection  coefficient  typical  of  turbulent  forced  convection  of  air  at 
moderate  velocity  over  a  flat  plate,  h  *  5,5  *  10-4  Btu/(ft2-s-R).  This  value 
should  only  be  regarded  as  a  rough  average  because  of  the  well  known  variabi¬ 
lity  of  h  with  many  factors  (surface  orientation,  flow  regime,  fluid  proper¬ 
ties,  etc.)  The  heat  transfer  between  a  cabin  surface  at  temperature  Ts  and  a 
gas  zone  at  temperature  Tg  is  then 

^cvn  *  "VTs  *  V  <4-19> 
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where  As  is  the  area  of  contact  between  the  surface  and  the  gas  zone.  The 
model  computes  the  fraction  of  each  lining  surface  In  contact  with  each  zone. 
When  a  surface  is  in  contact  with  both  the  upper  and  lower  zones  two  surface 
temperatures  are  used,  one  for  each  section  of  contact.  These  temperatureas 
are  Independent  since  no  lateral  heat  conduction  within  the  surface  material 
is  assumed.  The  method  of  computing  surface  temperature  is  presented  in 
Section  4.8. 


4.5.2  Radiation 

The  radiation  budget  for  the  gas  zone  is  the  difference  between 
the  absorbed  and  emitted  radiation.  To  compute  the  radiation  emitted  by  a 
zone  the  following  assumptions  are  made: 

•  Each  zone  emits  and  absorbs  radiation  as  a  gray  gas  with 
and  absorption  coefficient  which  is  proportional  to  the 
smoke  concentration. 

•  A  mean  beam  length  treatment  is  used  to  compute  the 
emissivity  of  each  zone. 

These  approximations  were  adoped  by  Quintiere  for  a  two-zone  model  of  room 
fires  (Reference  18).  The  emissivity  of  a  zone  is 

e  *  1  -  exp[-(osx^mj(p  +  (4-20) 

where  as  is  the  proportionality  factor  which  when  multiplied  by  the  zone  smoke 
concentration,  x’^,  and  density,  p,  yields  an  absorption  coefficient  repre- 

sentltl ve  of  the  absorption  and  emission  by  the  smoke.  In  the  units  of  smoke 
opacity,  employed  by  DACFIR3,  as  has  the  value  0.1054  ft2/particle.  The 
absorption  coefficient  kg  Is  an  added  factor  to  correct  for  gas  band  radiation 
as  suggested  by  Quintiere.  The  suggested  value  of  0.1  ft-1  is  used  for  kg. 

The  mean  beam  length,  L,  is  estimated  by  the  standard  method 
for  gas  volumes  of  a  rectangular  shape. 

L  *  3.6  (Zone  Volume/Asurf )  (4-21) 

where  Asurf  Is  the  surface  area  of  the  zone  volume. 

Total  radiation  emitted  by  a  gas  zone  Is  given  by 

*5 rout  *  ^surf^g*  (4-22) 

where  Tg  Is  the  zone  gas  temperature  and  a  Is  the  Stefan-Boltzmann  constant 

(4.761  x  10~13  Btu/ft2-sec-R. )  The  zone  also  absorbs  radiation  from  the  cabin 
surfaces  In  contact  with  it  and  from  the  other  zone  in  the  compartment  through 
the  thermal  discontinuity  plane.  This  absorbed  radiation  is  computed  by 

<Irin  *  e°(surfIcesAs,nTs  +  ]  {4'23) 

where  the  summation  is  over  all  cabin  lining  surfaces  in  contact  with  the  zone 
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gas,  is  the  area  of  contact,  and  Ts  Is  the  materials  surface  temperature. 
The  radiation  from  the  neighboring  zone  is  a  function  of  the  neighboring  zone's 
emissivity,  e..,  and  temperature  Tnz,  and  the  area  of  the  thermal  discontinuity 
plane,  A^.  Note  that  the  cabin  surface  materials  are  assumed  to  have  an 
emissivity  of  unity. 

4.6  GAS  LAW 

The  pressure,  temperature,  density,  and  composition  of  the  gas  in  each 
zone  are  related  in  DACFIR3  by  the  ideal  gas  law.  An  Important  simplifying 
approximation  in  the  modeling  of  enclosed  fires  by  the  zone  method  comes  from 
the  observation  that  the  hydrostatic  variation  in  pressure  between  the  floor 
and  ceiling  of  compartment  is  small  compared  to  the  absolute  pressure.  This 
Is  true  even  when  the  compartment  contains  a  layer  of  very  hot  gas.  Because 
the  change  In  pressure  with  height  is  small  compared  to  the  total  pressure,  we 
may  ignore  this  variation  and  use  a  single  reference  pressure  value  in  the  gas 
law  for  both  zones.  In  DACFIR3  we  choose  the  reference  pressure  to  be  the 

at  the  compartment  floor,  Pf.  Under  the  assumption  that  all  compartment  floor 
levels  are  the  same,  the  reference  pressure  is  then  at  the  same  geometric 
height  in  all  compartments.  Further  discussion  of  this  approximation  for 
pressure  may  be  found  in  Quintiere  (Reference  13). 

By  assumption,  only  the  five  major  gas  species  affect  the  mean  molecular 
weight  of  the  zone  gas.  Thus,  the  gas  law  for  each  zone  in  the  ith  compart¬ 
ment  is  written 

■  »K '  j,  (4-2,) 


where  and  T^  are  the  zone  gas  density  and  temperature,  R  is  the  universal 

gas  constant  (1545.0  ft-lbf/lbmole-R),  xj  is  the  mass  fraction  of  the  jth 
major  gas,  and  wmj  the  molecular  weight  of  the  jth  gas  species. 

4.7  ZONE  VOLUMES  AND  THE  THERMAL  DISCONTINUITY 

The  position  of  the  thermal  discontinuity  is  determined  from  the  fact 
that  the  total  volume  of  each  compartment  must  remain  fixed.  Thus  the  sum  of 
the  upper  and  lower  zone  volumes  must  equal  a  constant, 

V*  *  V.!  +  V1  *  constant  .  (4-25) 

t  u  X 

Since  the  cabin  cross  sectional  area  Is  fixed,  the  thermal  discontinuity 
height,  z^.  Is  a  function  of  only 

(4-26) 
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The  form  of  the  function  f  depends  upon  the  particular  shape  of  the  cabin 
cross  section.  For  example,  if  the  cross  section  is  a  simple  rectangle  of 
width  Cw  and  height  Ch  (4-26)  becomes 

z\  .  V«/(C’CW)  (4-27) 

where  is  the  compartment  length.  If,  on  the  other  hand,  th£  cross  section 

is  not  rectangular  but  contains  “steps"  as  in  Figure  2,  the  function  may  not 
have  a  convenient  closed-form  representation.  In  this  case,  an  implicit  eva¬ 
luation  of  (4-26)  is  employed  wherein  an  estimate  for  zi  is  used  with  the 
cabin  cross  section  profile  to  compute  a  trial  value  of  V^.  The  trial  value 
is  compared  to  the  current  true  value  which  is  known  from  (4-25).  An  itera¬ 
tive  correction  scheme  is  used  to  adjust  the  estimate  of  z^  until  it  results 
in  a  computed  V*  with  0.1  percent  of  the  true  value.  Further  discussion 

procedure  for  solving  all  of  the  cabin  atmosphere  equations  is  given  in 
Section  4.9. 


4.8  MATERIALS  SURFACE  TEMPERATURE 


Hot  gases  and  radiation  will  heat  the  surfaces  of  the  cabin  materials, 
especially  those  in  contact  with  the  upper  zone.  The  upper  heated  surfaces 
then  cut  down  the  heat  losses  from  the  upper  zone  gas  while  those  in  contact 
with  the  lower  zone  will  heat  the  lower  zone  gas.  DACFIR3  computes  the  sur¬ 
face  temperatures  of  the  cabin  lining  materials:  sidewalls,  ceiling  panels, 
floor  coverings,  etc.  These  structures  generally  consist  of  one  or  more 
lightweight  layers  of  polymeric  material  backed  by  thermal  insulation.  Since 
the  layers  are  thin,  a  lumped  thermal  model  should  be  an  adequate  approxima¬ 
tion  for  their  temperature  behavior. 


Each  part  of  the  cabin  lining  Is  considered  to  be  a  sheet  of  thickness 
as,  density  pm,  and  heat  capacity  c,,,.  The  total  surface  area  of  each  of  these 
sheets  is  fixed  an  the  fraction  of  that  surface  area  in  contact  with  each  gas 
zone  is  computed  at  every  integration  step  in  the  program.  The  important  heat 
transfer  mechanisms  are  taken  to  be  convection  to  the  zone  gas,  incoming  and 
emitted  radiation  at  the  Inside  (cabin  interior)  surface  and  conduction 
through  the  sheet  rear  face  in  contact  with  the  Insulation  layer. 


The  first  law  applied  to  a  unit  area  of  a  cabin  lining  material  yields 

V.“  £  ■  Vln  *  W  *  "<yTs>  *  T?n  <VT->  <4-28> 


where  Ts  is  the  materials  surface  temperature  (assumed  to  be  uniform 
throughout  the  sheet  thickness  as),  q^n  is  the  incoming  thermal  radiation, 
9rout  the  emitted  radiation  by  the  surface,  Tg  is  the  gas  temperature, 
k^n  is  the  insualtion  thermal  conductivity,  and  sin  is  the  insulation 


34 


thickness.  The  insulation  is  assumed  to  be  a  simple  linear  thermal  resistance 
of  negligible  heat  capacity.  Its  outside  is  assumed  to  always  be  at  the 
ambient  temperature,  I..  The  incoming  thermal  radiation  is  that  emitted  by 
both  gas  zones,  as  discussed  in  Section  4.5.2,  while  that  emitted  by  the  sur¬ 
face  is  qp0ut  -  oT$4  .  The  values  of  k-jn  and  s^n  are  supplied  by  th*  user  for 
each  lining  service. 


4.9  NUMERICAL  SOLUTION  OF  THE  ATMOSPHERE  EQUATIONS 


Equations  (4-1)  and  (4-2)  give  the  rate  of  change  of  the  mass  of  gas  in 

each  zone  of  each  compartment,  and  M1.  respectively.  The  amounts  of  mass 

u  » 

are  not  as  convenient  to  work  with  as  are  the  zone  densities  given  by 


•>1 =  MX 


1  u1  /t/1 

p  *  M  /V 
Hu  u  '  u 


(4-29) 

(4-30) 


With  these  relationships  the  model  of  the  cabin  atmosphere  is  complete.  Table 
2  summarizes  the  dependent  variables  of  the  model  and  the  primary  equations 
which  relate  them.  Briefly,  the  equations  employed  are  the  conservation  rela¬ 
tionships  for  mass,  species,  energy,  and  total  compartment  volume;  the  gas  law 
for  both  zones;  conservation  of  energy  for  the  materials;  the  definition  of 
density;  and  the  geometric  relationship  of  the  boundary  between  the  zones  to 
the  lower  zone  volume  and  cabin  cross  section.  All  the  symbols  in  Table  2 
carry  the  superscript  i  to  indicate  values  for  the  ith  compartment.  When  more 
than  one  compartment  is  present  there  is  a  set  of  (8+j+k)  variables  and 
(8+j+k)  equations  for  each  compartment.  The  entire  set  of  i  •  (8+j+k) 
equations  is  coupled  by  the  flows  through  the  vents  connecting  compartments. 


The  numerical  solution  of  this  large  set  of  mixed  algebraic  and  ordinary 
differential  equations  poses  a  formidable  problem.  The  earlier  versions  of 
DACFIR  useu  simple  one-step  Euler  integration  for  the  differential  equations 
with  the  necessary  algebraic  relationships  applied  explicitly  where 
appropriate.  While  this  method  was  simple  and  straightforward,  the  time  steps 
required  for  a  stable  solution  were  very  small,  on  the  order  of  0.001  second. 
The  non-linear  character  of  the  equations,  particularly  the  set  used  in 
Version  2  to  compute  flow  through  vents,  made  the  numerical  solution  margi¬ 
nally  stable  at  even  the  smallest  step  sizes. 

A  superior  method  for  this  problem  Is  the  implicit  trapezoidal  rule 
Integration  used  by  Emmons,  Mltler,  and  Trefethen  (Reference  19)  for  the 
Harvard  Computer  Fire  Code.  This  method  is  well  explained  in  the  reference  so 
that  only  a  brief  summary  need  be  given  here. 

Each  of  the  differential  equations  of  the  governing  set  can  be  written 
as 


d 

dtxn 


f((x}) 


(4-31) 
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TABLE  2 


SUMMARY  OF  VARIABLES  AND  EQUATIONS  OF  THE  CABIN  ATMOSPHERE  MODEL 


Variable 

Symbol 

Obtained  from  Equation(s) 

Lower  zone  species 
mass  fractions 
(j  values)* 

i 

x£j 

(4-3)  (j  equations) 
and  (4-1) 

Upper  zone  species 
mass  fractions 
(j  values) 

i 

xuj 

(4-4)  (j  equations) 
and  (4-2) 

Pressure 

(4-24) 

Lower  zone  density 

i 

pji 

(4-29)  and  (4-1) 

Upper  zone  density 

•1 

(4-30)  and  (4-2) 

Lower  zone  temperature 

(4-8)  and  (4-2) 

Upper  zone  Temperature 

(4-9)  and  (4-1) 

Lower  zone  volume 

*1 

(4-25) 

Upper  zone  volume 

V1 

u 

(4-24)  and  (4-30) 

Thermal  discontinuity 

(4-26) 

position 

Materials  surface 
temperature 

4 

(4-28)  (k  equations)** 

*  Minimum  value  of  j  is  5  and  the  maximum  11 

'*  Minimum  value  of  k  is  6  and  the  maximum  22  per  compartment 
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where  x»  is  one  of  the  set  {x}  of  dependent  variables.  In  the  trapazoidal 
rule  integration,  the  value  of  xn  at  time  t+At  is  estimated  by  the  formula 

xn(t+At)  «  xn(t)  +  l/2[f (t)  +  f (t+At)]  (4-32) 

Equation  (4-32)  is  an  Implicit  relation  for  xn(t+At)  since  f(t+At)  involves 
the  specification  of  the  set  (x)  at  t+At.  When  the  equations  represented  by 
(4-32)  are  combined  with  the  original  algebraic  equations  the  problem  of  pre¬ 
dicting  the  values  of  the  variables  at  t+At  from  the  values  at  t  becomes  that 
of  solving  a  large  set  of  non-linear  algebraic  equations. 

A  number  of  methods  are  available  to  solve  non-linear  algebraic 
equations,  each  with  its  particular  advantages  and  drawbacks.  For  DACFIR3  we 
have  chosen  the  Newton-Raphson  method  because  of  its  general  reliability.  The 
Newton-Raphson  method  is  an  iterative  technique  in  which  an  initial  guess  for 
the  solution  {x>  at  t+At  is  repeatedly  corrected  using  the  equation  set  {f}. 
Each  new  correction  (iteration)  produces  a  new  estimate,  tx}^,  of  the  solu¬ 
tion.  When  the  (k+1)  solution  is  sufficiently  close  to  the  previous  solution, 
k,  the  set  is  considered  solved.  The  closeness  of  two  solution  sets  is  judged 
by  testing  the  difference  between  the  corresponding  values  in  the  two  sets 
against  a  pre-set  tolerance. 

The  algorithm  of  the  Newton-Raphson  method  involves  writing  Equation 
(4-32)  as 


{x(t+At) }  =  (g({x(t+At)>,  (x(t)})},  (4-33) 

where  the  brace  notation  is  used  to  indicate  the  complete  system  of  n 
equations.  The  set  of  n  functions  g  involve  the  variables  at  t  and  t+At.  A 
simple  rearrangement  of  (4-33)  and  the  definition  of  another  function,  F,  as 
Fn  s  9„  -  ^n (t+At )  gives 

{F ( (x(t+At) } ,  {x(t)})>  -  0  .  (4-34) 

Each  of  the  n  functions  F  is  Identically  zero  when  x(t+At)  solves  the  equation 
set.  Since  the  values  of  x  at  time  t  are  known  and  do  not  change  during  the 
iteration,  we  may  regard  the  set  {F}  as  functions  only  of  the  estimates  of  x 
at  t+At,  that  is,  during  one  time  step  we  wish  to  solve  for  {x}  such  that 

(F(  {x} )}  ■  0  (4-35) 

where  it  Is  now  understood  that  tx)  Is  evaluated  at  t+At. 

The  equation  for  generating  each  successive  estimate  of  the  solution  is 

{x}k+l  *  {x}k  -  C03”1{F({x>k)>  (4-36) 

where  [J]  is  the  n-by-n  Jacobian  matrix  of  the  set  F.  Each  element  of  the 
Jacobian  matrix  Is  defined  as 

3F1 
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Following  the  method  used  by  Emmons,  et  al.,  the  elements  of  the 
Jacobian  are  evaluated  numerically  from  the  definition  of  the  partial 
derivative. 


'U 


My")  -W 


(4-38) 


where  h  is  a  small  quantity  added  to  the  jth  variable  of  the  set  {x}  while  the 
other  members  of  the  set  are  held  fixed. 

Finally  to  solve  (4-36)  for  {xj^+i,  we  can  rearrange  once  more  to  obtain 


[J]({x}k  -  {x}k+1)  =  (F>  .  (4-39) 

This  way  of  writing  the  correction  equation  allows  us  the  use  of  standard 
pre-tested  computer  routines  for  solving  systems  of  linear  equations. 

The  numerical  technique  described  above  is  best  accomplished  when  all 
the  variables  are  of  about  the  same  order  of  magnitude,  say  on  the  order  of 
one.  In  DACFIR3  a  scaling  procedure  Is  adopted  from  (Reference  19)  which  sca¬ 
les  the  variables  to  0(1)  before  the  evaluation  of  the  Jacobian  and  the  solu¬ 
tion  of  (4-39).  The  new  guess  at  the  solution  can  then  be  compared  to  the 
previous  guess,  also  given  in  scaled  magnitude,  using  the  same  tolerance  for 
the  difference  between  each  pair  of  variables.  If  further  iteration  is  needed 
the  variables  are  rescaled  back  to  their  physical  magnitudes  for  the  eva¬ 
luation  of  the  functions  {F}. 

One  last  point  concerning  the  numerical  solution  is  that  the  variables 
and  equations  for  the  cabin  lining  surface  temperatures  are  not  included  in 
the  iterative  procedure.  The  thermal  inertia  of  the  materials  is  large  enough 
to  keep  the  derivatives  on  the  left  in  (4-25)  small  during  a  time  approxima¬ 
tion  to  keep  the  surface  temperatures  constant  during  the  iterative  time  step 
Is  reached,  update  the  surface  temperatures  by  single  step  Euler  Integration. 
This  approximation  removes  from  six  to  20  variables  and  equations  for 
every  compartment— obviously  a  significant  simplification. 
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SECTION  5 
RADIATION 


» 


Thermal  radiation  to  the  fuel  surface  is  the  single  governing  parameter 
for  determining  ignition,  flame  spreading,  product  emission  rates,  and  dura¬ 
tion  of  burning  for  the  cabin  materials.  The  total  incoming  radiation  at  an 
element's  surface  is  taken  to  consist  of  two  parts:  that  Ijrom  nearby  flames, 
and  that  from  the  hot  combustion  products  which  form  the  upper  gas  zone.  We 
term  the  first  part  the  local  flame  radiation  and  the  second  part  the  gas  zone 
radiation.  Local  flame  radiation  is  computed  by  making  a  number  of 
simplifying  assumptions  about  flame  structure,  radiation  temperature,  and 
optical  properties.  This  development  is  given  in  Section  5.1  below.  Gas  zone 
radiation  has  been  discussed  above  in  Section  4.5  where  the  net  radiation 
balances  for  the  gas  control  volume  were  given.  In  Section  5.2  we  show  how 
the  zone  radiation  to  a  material  element  is  estimated  from  the  zone  emittance 
and  the  element  position. 


5.1  FLAME  RADIATION 

Three  estimates  of  the  radiation  from  the  flames  of  a  fire  are 
made  by  DACFIR3,  each  for  a  separate  purpose.  The  flux  to  elements 
just  outside  of,  but  adjacent  to,  the  fire  base  area  is  qj,  the 

average  flux  to  elements  forming  the  fire  base  area  is  q2,  and  the 

flux  to  elements  outside  but  in  the  vicinity  of  the  tip  of  the  flames 

2 

is  q3.  All  these  fluxes  are  given  in  units  of  Btu/ft  •  sec.  Figure 
10  illustrates  these  three  quantities. 

////// 1  ^/// 
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To  compute  the  radiation  fluxes,  the  following  assumptions  are  made 
about  the  flames: 

•  The  flame  volume  may  be  approximated  as  a  right  circular  cylinder  of 
height  hf  and  base  diameter  2y0. 

•  The  flame  volume  radiates  as  a  gray  gas  with  a  spatially  uniform  tem¬ 
perature  and  absorption  coefficient. 

For  sooty  flames,  the  type  expected  In  cabin  fires,  practically  all  of 
the  radiation  Is  due  to  the  soot  pa  tides  In  the  flame.  Since  measurements 
of  the  smoke  generated  by  the  cabin  materials  are  available,  the  emlsslvitles 
of  fires  of  these  materials  can  be  estimated  from  the  smoke  generation  data. 
Under  a  number  of  simplifying  assumptions  a  method  was  developed  to  express 
the  absorption  coefficient  in  terms  of  the  smoke  generation  rate  as 

kf  •  0. 21  p H  (hfgf 1/2  (5-1) 

where  k#  is  the  flame  absorption  coefficient  (ft-1),  and  pM  Is  the  smoke 
generation  rate  of  the  fuel  (particles/ft2  .  sec).  The  derivation  of  equation 
(5-1)  Is  given  in  Appendix  0  of  Reference  2.  Using  this  estimate  of  the 
flame  absorption  coefficient,  the  emlttance  at  the  fire  base  center  Is 

oc  -  [1  -  exp(-1.8kfhf)3  +  [1  -  exp(-1.8kfy0)] 

(5-2) 

-  {1  -  exp[-1.8kf(hf2  +  y02)l/2]>. 

where  y0  Is  the  flame  cylinder  base  radius  (ft).  This  expression  was  devel¬ 
oped  by  Dayan  and  Tien  (Reference  20)  for  a  cylindrical  flame. 

Dayan  and  Tien's  analysis  shows  that  the  radiation  flux  at  the  edge  of 
the  flame  base  Is  about  one  half  of  the  value  at  the  base  center  for  flames 
with  a  helght-to-diameter  ratio  on  the  order  of  two.  Since  this  Is  a  typical 
value  for  the  flame  lengths  of  cabin  materials  we  express  as 

qi  »  0.50^,  .  (5-3) 

The  flame  emissive  power,  et>.  Is  taken  to  be  a  constant  for  all  fires.  A 
study  by  Fu  (Reference  21)  of  the  radiation  from  aircraft  fuel  fires  of  base 
diameters  of  several  feet  Indicates  that  the  effective  emissive  power  for 
these  fires  Is  about  16.3  Btu/ft2  •  sec.  A  review  of  data  on  smaller  scale 
fires  (about  1  ft  base  diameter)  of  plastic  materials  by  deRis  (Reference  10) 
shows  that  the  emissive  power  of  these  flames  Is  In  the  range  of  10  to  25 
Btu/ft2  •  sec.  In  view  of  this,  the  value  16.3  Is  used  for  e^  In  DACFIR3. 

For  q?»  the  average  flux  over  the  base,  the  analysis  of  Dayan  and  Tien 
Indicates  that  the  average  vaue  Is  approximately  84  percent  of  the  value  at 
the  center.  Therefore  we  write  q£  as 

q2  ■  0.84oceb  .  (5-4) 

The  heat  flux  to  an  element  on  a  surface  above  a  fire  Is  estimated  as 


40 


(5-5) 


q3  -  1/2[1  -  (4C2  -  3 )//  (4c2  +  9)(4C2  +  l)]eb  . 

This  expression  Is  obtained  by  using  the  circular  cylinder  model  of  the  flame 
and  taking  the  target  element  to  be  located  a  distance  z  above  the  tip  and 
one-half  the  radius  (y0)  from  the  axis  of  the  cylinder.  The  ration  z/y0  Is 
the  factor  c  of  (5-5). 

The  smoldering  range  xp,  Is  the  distance  away  from  the  flame  foot  at 
which  the  radiation  drops  toKthe  level  qp.  To  compute  xp  we  again,  employ  the 
analysis  of  Dayan  and  Tien,  first  expressing  xp  as  v 

*p  ■  x  -  y0 

where  x  Is  the  distance  from  the  flame  cylinder  axis  to  the  point  where  q  * 
qp.  The  value  of  x  Is  obtained  by  solving  the  cubic  equation 

wyQx3  +  (0.5hf 2  -  »y02)x2  -  u02hf2(eb/qp  -  0.5)  .  (5-7) 

Of  the  three  roots  of  (5-7),  the  largest  real  root  Is  always  the  physically 
correct  value.  The  computer  code  Is  designed  to  select  this  value  In  all 
cases. 

The  flame  height,  hf.  Is  an  Important  parameter  In  all  of  the 
expressions  above.  In  DACFIR3  the  flame  height  calculation  method  of 
Fang  (Reference  15)  Is  adopted  with  a  slight  modification.  Fang's  analysis  of 
freely  burning  fires  has  shown  that  the  visible  flame  height  correlates  reso- 
nably  well  with  the  point  on  the  plume  axis  where  the  temperature  of  the  com¬ 
bustion  products  has  fallen  to  about  400°C  (752°F)  or  about  2.25  times  the 
room  temperature.  His  relationship  for  the  flame  height  Is 

h,  -  (1.49  +  0.916  K*/5)Pa1/5N  2/5y„  (5-8) 

T  d  a  D  0 


where 

Ka  3  <Ec/Ep)4(1  ‘  “)C2*25fl  +  “K’V  "  1)/P0,]3/Cl*95“3f23{1  '  ps‘)L 
Pa  *  w^2^Ep4^1  " 

Nb  3  P0uo/(p/^o)* 

fl  *  -  po')/po'  +  y,\’ 

f2  -  m/Po1  +  Y/xi  . 

^  -4 

po'  3  p0/p*  To'  3  T0/T  *  an<^  Ps*  "  ps/p- 

It  can  be  seen  that  this  formidable  expression  Involves  porpertles  of  the 
fuel:  p„,  u0,  y,  w,  T0;  of  the  exterior  (to  the  fire)  gas,  T1,  p;  and  of 
the  base  area,  y0.  In  addition,  the  empirical  entrainment  constants  Ec  and 
Ep  appear.  Comparing  (5-8)  to  experiments  on  small  scale  fires  of  n-heptane 
Fang  found  good  agreement  using  Ec  -  0.25  and  Ep  -  0.05.  When  compared  to  the 
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data  on  plastic  fires  of  Modak  and  Croce  (Reference  22),  we  found  that  the 
values  of  0.3  for  Ec  and  0.25  for  Ep  produced  better  predictions  of  the 
observed  flames.  We  have  thereof re  adopted  these  larger  values  for  the  esti¬ 
mation  of  flame  lengths  in  DACFIR3.  When  entrainment  is  calculated  for  the 
mass  conservation  equations  of  the  cabin  atmosphere  model  the  smaller  values 
for  Ec  and  Ep  given  in  Section  4.3  are  retained. 

5.2  GAS  ZONE  RADIATION 

In  principle,  a  material  element  receives  incoming  thermal  radiation 
from  both  gas  zones  and  all  other  material  elements  that  can  be  viewed  from 
the  surface  of  the  element.  However,  conditions  during  the  early  stages  of  a 
cabin  fire  are  such  that  only  the  radiation  from  the  hot  upper  zone  is 
significant.  DACFIR3  includes  two  calculations  of  the  radiation  flux  to 
material  elements  from  the  upper  zone  gas:  radiation  to  elements  which  are  in 
direct  contact  with  the  upper  zone,  and  radiation  which  travels  through  the 
lower  zone  before  reaching  elements  in  contact  with  the  lower  zone  gas.  Both 
of  these  radiation  fluxes  are  added  to  the  flux  levels  described  in  the  pre¬ 
vious  section  to  determine  the  material's  fire  behavior. 

The  computation  for  materials  in  the  upper  zone  uses  the  same  gray  gas 
and  mean  beam  length  approximations  as  described  in  Section  4.5.2.  The  energy 
arriving  per  unit  area  per  unit  time  is 

qzu  a  [1  -  exp(-kuLu)]oTu4  (5-9) 

where  ku  is  the  upper  zone  absorption  coefficient  and  Lg  is  the  mean  beam 
length.  Since  the  element  surface  Is  in  contact  with  tne  zone  gas  the  view 
factor  is  unity. 

For  lower  zone  elements  the  radiation  intensity  from  the  upper  zone  may 
be  reduced  by  absorption  in  the  lower  zone  and  by  the  view  factor  between  the 
element  and  the  lower  surface  of  the  upper  zone,  i.e. ,  the  thermal  discon¬ 
tinuity  plane.  The  flux  Is  written 

qzt  *  exP(-k*Li)FCl  -  exP(-kuLu)]oTu4  (5-10) 

where  k»  and  LA  are  the  absorption  coefficient  and  beam  length  for  the  lower 
zone.  The  value  of  the  view  factor  F  depends  on  the  position  of  the  thermal 
discontinuity  and  the  location  of  the  element.  To  compute  a  separate  value  of 
F  for  each  lower  zone  element  Is  Impractical  for  the  number  of  elements 
involved  so  a  representative  view  factor  is  used.  This  value  is  that  for  an 
element  facing  upward  at  the  center  of  the  cabin  floor, 

F  *  (2/ir) [ (a/A)tan“ 1  ( b/A)  +  (b/BJtan’1  (a/B)]  (5-11) 

where  a  *  Cw/(  zd),  A  -  (1  +  a2)1/2,  b  =  CL/(2zd),  B  *  (1  +  b2)1/2, 

Cy  is  the  compartment  width  at  the  floor,  C[_  is  the  compartment  length,  and 
zd  is  the  thermal  discontinuity  height. 
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SECTION  6 

SAMPLE  RESULTS  OF  THE  SIMULATION  OF  SEVERAL  CABIN  FIRE  TESTS 


This  section  presents  selected  results  from  the  simulation  by  DACFIR3  of 
three  full-scale  fire  tests.  These  results  are  presented  to  demonstrate  the 
cpability  of  the  model  to  reproduce  test  results,  given  the  best  available 
input  data.  We  have  selected  two  of  the  most  Important  output  quantities  of 
DACFIR3,  the  predictions  of  gas  temperature  and  species  concentration,  to  com¬ 
pare  to  the  experimental  results.  Many  other  quantities  are  Included  In  the 
model's  output,  a  complete  list  of  which  is  given  in  Volume  II. 

The  cabin  fire  tests  discussed  here  are  selected  from  a  large  series  of 

tests  conducted  by  Bricker,  et  al.  (Reference  23)  in  a  B- 737  fuselage  spe¬ 

cially  prepared  and  instrumented  for  repeated  full  scale  testing.  The  cabin 
cross  section  measured  approximately  11  feet  wide  at  the  floor  and  seven  feet 
from  floor  to  ceiling  at  the  highest  point.  Moveable  bulkheads  (partitions) 
allowed  for  testing  at  two  cabin  lengths:  56  feet  fore-to-aft  giving  a  total 

cabin  volume  of  about  3700  cubic  feet;  and  20  feet  in  length  with  a  total 

volume  of  1320  cubic  feet.  The  interior  lining  of  the  cabin  was  lined  with 
steel  and  aluminum  sheet  backed  by  high  temperature  resistant  insulation.  The 
lining  configuration  simulated  a  modern  standard-body  cabin  using  overhead 
pull-down  stow  bins. 

Instrumentation  consisted  of  an  array  of  thermocouple  trees  for 
measuring  air  temperature,  several  gas  sampling  ports  and  automatic  gas  analy¬ 
sis  equipment,  smoke  meters,  flowmeters,  heat  flux  gages,  and  load  cells  for 
measuring  rate  of  weight  loss  during  burning.  Ventilation  of  the  cabin 
occurred  through  doors  In  the  fore  and  aft  bulkheads.  For  most  tests  a  large 
blower  supplied  a  constant,  known  air  flow  directed  Into  the  cabin  through  the 
forward  bulkhead  door.  In  order  to  simulate  a  wind  from  an  open  door.  Several 
tests  were  conducted  without  forced  ventilation,  the  bulkhead  doors  being  open 
to  the  exterior.  This  condition  Is  termed  natural  ventilation. 

The  three  tests  simulated  by  DACFIR3  are  summarized  below. 

Test  3B.  In  this  test  a  pan  of  aircraft  turbine  fuel  (Jet  A)  was  burned 
in  the  bare  (unfurnished)  cabin  of  the  56  ft  length.  The  objective  of  this 
test  was  to  study  the  behavior  of  pan  fires  of  fuel  which  were  used  as  igni¬ 
tion  sources  in  later  tests  Involving  materials.  Located  at  floor  level  In 
the  center  of  the  cabin,  the  square  pan,  four  square  feet  in  area,  held  4.5 
liters  of  fuel.  A  load  cell  measured  the  weight  loss  rate  of  the  pan  during 
the  burning  period,  which  lasted  about  10  minutes.  Natural  ventilation  was 
used  in  this  test.  The  dimensions  of  the  bulkhead  doors  were  five  feet  high 
by  two  and  one  half  feet  in  width.  The  tops  of  the  doors  were  five  feet  nine 
inches  above  the  floor,  leaving  small  sills  nine  Inches  high  at  the  floor. 
Seven  thermocouple  trees  spaced  at  roughly  equal  intervals  fore-to-aft 
measured  air  temperature  continuously  at  20,  40,  60,  and  75  Inches  above  the 
floor.  Measurements  of  gas  composition  were  made  at  several  horizontal  posi¬ 
tions  30  and  60  inches  above  the  floor  at  45  second  Intervals. 

Test  5A.  Conditions  In  this  test  were  Identical  to  those  in  Test  3B 
with  two  exceptions.  First,  the  cabin  length  was  20  feet  Instead  of  56  feet. 
Second,  forced  ventilation  was  employed.  An  airflow  of  500  cubic  feet  per 
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minute  (cfm)  entered  through  the  forward  bulkhead  door,  the  aft  door  being 
open  to  the  outside.  Because  of  the  shorter  length  of  the  cabin,  only  five 
equally  spaced  thermocouple  trees  were  used  for  air  temperature.  The  vertical 
positions  of  the  thermocouples  and  gas  sampling  ports  were  the  same  as  In  3B. 

Test  14A.  This  test,  conducted  using  the  56-foot  long  cabin,  involved  a 
seat  row  and  celling,  sidewall,  and  passenger  service  unit  (PSU)  materials.  A 
one-square  foot  pan  of  Jet  A  fuel  burned  beneath  the  outboard  seat  of  a  row  of 
three  surplus  aircraft  seats.  Seat  construction  was  of  a  typical  wool -nylon 
blend  fabric  covering  fire-retarded  polyurethane  foam  padding.  On  the  wall 
next  to  the  outboard  seat  a  panel  of  Tedlar/epoxy-flberglass/Nomex  honeycomb 
sidewall  material  extended  from  the  floor  to  the  PSU  surface.  The  PSU  surface 
was  simulated  by  a  sheet  of  polycarbonate  stock.  Both  the  sidewall  panel  and 
the  PSU  sheet  were  four  feet  In  length  (fore-aft  direction).  Celling  panel 
material  of  the  same  construction  as  the  sidewall  covered  a  four  foot  length 
on  the  stow  bln  face  surface  and  a  20  foot  length  on  the  celling  proper.  The 
bulkhead  doors  had  the  same  dimensions  as  In  Tests  38  and  5A.  The  ventilation 
was  500  cfm  through  the  forward  door. 

For  simulation  of  these  tests  by  DACFIR,  Input  data  was  collected  from 
a  number  of  sources.  Some  of  these  quantities,  such  as  the  cabin  dimensions, 
ventilation,  and  materials  location,  were  taken  directly  from  the  test 
description.  One  quantity,  the  mass  burning  rate  of  the  ignition  (fuel  pan) 
fire,  was  obtained  from  the  measured  weight  loss  rates  of  the  pan  In  each 
test.  The  materials  flammability  data  for  Test  14A  was  selected  from  the  data 
set  collected  for  the  DACFIR2  validation  program.  This  data  set  is  described 
.  In  (Reference  2).  Although  flammability  data  was  not  obtained  from  materials 
Identical  to  those  used  In  the  tests,  the  data  was  obtained  from  materials 
that  were  nominally  the  same.  In  terms  of  chemical  composition  and  physical 
construction.  We  believe  that  the  materials  data  used  Is  sufficiently  charac¬ 
teristic  of  the  materials  used  In  the  full-scale  tests  to  make  the  comparison 
meaningful.  Finally,  a  number  of  Input  quantities,  such  as  the  heat  of  com¬ 
bustion  of  the  Ignition  source  fuel  and  the  thermal  properties  of  the 
materials,  were  taken  from  the  general  literature  or  given  "best  estimate" 
values.  In  some  cases  the  accuracy  of  these  estimates  Is  very  good.  In  other 
cases,  however,  there  Is  considerable  question  as  to  the  appropriate  values. 
When  judging  the  performance  of  the  model  In  comparison  to  test  results,  one 
must  keep  In  mind  the  problem  of  providing  the  best  estimates  of  the  necessary 
Input  data. 

Figure  11  shows  the  comparison  of  the  predicted  and  measured  gas  tem¬ 
perature  for  Test  38.  The  experimental  curve  was  obtained  by  a  simple  average 
of  the  readings  of  the  seven  thermocouples  at  60  Inches  and  the  seven  at  75 
Inches.  Thermocouple  readings  were  digitized  at  one  second  Intervals  although 
only  the  values  at  each  30  seconds  have  been  plotted  In  the  figure.  The  pre¬ 
dicted  curve  Is  the  upper  zone  temperature  from  DACFIR3  computed  using  the 
measured  mass  loss  rate  of  the  fuel  pan  In  the  test.  This  rate  was  a  very 
constant  0.00311  lbm/ft^-sec  throughout  the  first  six  minutes  of  the  burn. 

The  rate  was  converted  to  a  rate  of  heat  release  by  assuming  an  effective  heat 
of  combustion  for  Jet  A  fuel  of  16,650  Btu/lbm.  The  effects  of  combustion 
Inefficiency  and  sensible  heat  loss  by  flame  radiation  are  taken  Into  account 
through  this  effective  value. 
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Figure  11.  Predicted  and  Measured  Upper  Zone  Gas  Temperature  for  Test  3B. 

The  figure  shows  fair  agreement  between  the  prediction  and  the  test,  the 
maximum  difference  being  about  80°F  at  300  seconds.  The  intitial  rate  of  rise 
of  the  DACFIR  curve  is  partially  determined  by  a  feature  of  the  model  which 
allows  the  user  to  "ramp-in"  the  effect  of  the  ignition  fire  over  a  specified 
period.  During  this  period,  the  rates  of  heat,  smoke,  and  gas  release  are 
increased  linearly  from  zero  to  the  constant,  maximum  values  given  in  the 
input.  From  observation  of  the  test  films,  a  ramp-in  time  of  20  seconds  was 
selected  for  all  simulations. 

After  the  ramp-in  period,  the  figure  shows  the  predicted  temperature 
dropping  slightly  during  the  first  two  minutes  and  then  rising  slowly  there¬ 
after.  This  behavior  reflects  the  influence  of  the  slow  heating  of  the  walls 
and  the  subsequent  decrease  of  heat  loss  from  the  upper  zone. 

The  upper  zone  thickness  computed  by  DACFIR  for  this  test  was  quickly 
established  at  about  2.5  feet  after  the  ramp-in  period.  It  remained  close  to 
this  value  throughout  the  simulation.  (This  value  of  the  zone  thickness  was 
used  to  select  which  test  thermocouples  would  be  averaged  for  the  comparison.) 

Another  comparison  of  predicted  and  measured  temperature  of  the  accumu¬ 
lating  hot  gases  is  shown  In  Figure  12.  Test  5A  differed  from  3B  in  that  the 
cabin  volume  was  only  about  one  third  of  that  in  3B  and  that  a  fixed  500  cfm 
forced  ventialtlon  was  used.  The  effect  of  the  smaller  vaolume  Is  readily 
apparent  in  the  figure  where  both  the  predicted  and  measured  values  reach 
about  500°F.  In  contrast  to  38,  the  prediction  of  temperature  In  5A  is  higher 
than  that  measured.  Again  a  measured  burning  rate  from  the  test,  0.0058 
lbm/ft^-sec,  was  used  as  Input  to  the  simulation. 

Figure  13  gives  the  predicted  and  measured  gas  temperature  in  Test  14A, 
the  test  involving  materials.  Figures  14  and  15  show  the  comparison  of  the 
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12.  Predicted  and  Measured  Upper  Zone  Gas  Temperature  for  Test  5A. 


predictions  and  measurements  of  three  toxic  gases,  CO,  HF  and  HCN.  All  three 
figures  show  fair  agreement  between  prediction  and  measurement  for  the  first 
four  or  five  minutes.  During  this  time  the  primary  heat  source  is  the  igni¬ 
tion  source  fire  burning  at  0.00259  1bm/ft2-sec.  The  pan  size  in  this  test 
was  one  square  foot.  As  burning  develops  on  the  seat  and  sidewall  the 
emissions  of  the  toxic  trace  gases  increase. 

The  model's  prediction  of  CO  accumulation  Is  primarily  ahead  of  the 
measurement  while  the  predictions  of  HF  and  HCN  are  mostly  behind.  The  noti¬ 
ceable  recession  in  the  predicted  gas  concentrations  starting  around  four 
miunutes  and  lasting  for  about  30  seconds  is  the  result  of  decrease  in  the 
total  emission  rates  during  this  period.  This  decrease  seems  to  be  a  con¬ 
sequence  of  the  way  in  which  the  model  selects  fire  base  areas  and  thus  flame 
volumes,  flame  radiation  levels,  and  subsequently,  emission  rates. 

The  situation  then  reverses  after  five  minutes  when  the  predicted  gas 
concentrations  for  CO  and  HCN  begin  to  exceed  the  measurement  and  the  HF  pre¬ 
diction  rises  rapidly,  closing  in  on  the  measured  value  at  six  minutes.  For 
temperature,  however,  the  model’s  prediction  is  well  below  the  measurement 
after  five  minutes.  The  origin  of  the  rapid  temperature  rise  In  the  test 
results  is  most  probably  an  acceleration  of  the  burning  of  the  seats.  The 
seat  fire  In  the  DACFIR  simulation  is  also  growing  steadily  during  this  time 
but  apparently  not  at  the  same  rate  as  in  the  test. 


CO  Concentration  (ppm) 


Figure  14.  Predicted  and  Measured  Upper  Zone  Carbon  Monoxide  Concentration 
for  Test  14A. 
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Figure  15.  Predicted  and  Measured  Upper  Zone  Hydrogen  Cyanide  and 
Hydrogen  Fluoride  Concentrations  for  Test  14A. 
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SECTION  7 
CONCLUSIONS 


Version  3  of  DACFIR  has  been  created  as  a  refinement  and  generalization 
of  the  cabin  fire  model.  The  major  differences  between  this  version  and  the 
earlier  versions  of  DACFIR  are  the  revised  cabin  atmosphere  model  which  allows 
for  multiple  compartments  and  the  prescribed  entry  of  exterior  fire  gases  and 
the  implicit  Newton-Raphson  integration  technique  for  solution  of  the  cabin 
atmosphere  equation  set.  Major  revisions  have  been  made  to  the  computer  code 
to  provide  a  more  modular  structure  making  the  code  easier  to  understand, 
maintain,  and  upgrade. 

Our  experience  in  using  DACFIR3  for  the  simulation  of  three  full-scale 
cabin  fire  tests  has  lead  to  the  following  conclusions: 

(1)  The  implicit  numerical  technique  for  solution  of  the  cabin 
atmosphere  equations  has  proven  to  be  superior  to  earlier  methods  in  terms  of 
both  stability  and  computer  time  required.  Some  initial  trial  and  error 
experimentation  is  required  to  find  the  optimum  combination  of  integration 
time  step,  convergence  tolerance,  and  maximum  number  of  Iterations  for  each 
case.  Once  this  combination  Is  Identified  the  model  appears  to  operate  with 
reasonable  efficiency  and  accuracy. 

(2)  The  cabin  atmosphere  model  performs  adequately  in  predicting  the 
average  upper  zone  temperature  for  fires  of  constant  rate  of  heat  release. 
Predictions  of  C02,  not  shown  in  this  report,  are  consistently  lower  than 
measured  in  the  test.  This  problem  is  probably  due  to  the  oversimplified  com¬ 
bustion  chemistry  assumed  by  the  mode?  which  reflects  the  current  lack  of 
understanding  of  the  fundamental  behavior  of  turbulent,  diffusion-controlled 
burning  of  liquids.  Smoke  and  trace  gas  accumulation  for  the  fuel  pan  fire 
were  not  compared  to  test  results. 

(3)  When  cabin  furnishing  materials  were  included  in  the  fire,  the 
model  performed  adequately  in  predicting  the  early  development  of  the  fire, 
both  for  temperature  and  trace  gas  accumulation.  As  the  fire  growth 
accelerates,  the  model's  predictions  generally  lag  the  test  results.  The 
apparent  explanation  for  this  disagreement  concerns  the  range  of  heat  flux 
levels  covered  by  the  materials  flammability  data.  This  data,  which  includes 
Ignition  delay  times,  flame  spread  rates,  and  rates  of  heat  and  gas  emission, 
was  measured  in  the  laboratory  at  imposed  heat  flux  levels  of  up  to  about  4.5 
Bcu/ft^-sec.  0ACFIR3's  computation  of  the  flux  levels  at  a  flame  base, 
however,  can,  and  often  are,  much  higher  than  this  value— sometimes  reaching 
12  or  13  Btu/ft‘-sec.  To  find  emission  rates  at  these  high  flux  levels,  the 
model  linearly  extrapolates  from  the  known  rates.  It  Is  quite  likely  that  the 
true  rates  of  emission,  etc.  at  these  high  fluxes  are  not  equal  to  the 
linearly  extrapolated  values.  In  particular,  we  might  expect  the  high  flux 
rates  of  heat  release  to  be  greater  than  the  extrapolated  values  while  some 
gas  emission  rates,  such  as  CO,  might  be  less.  Relatively  higher  heat  release 
rates  and  lower  gas  emission  rates  would  result  in  higher  upper  zone  tem¬ 
peratures  and  lower  trace  gas  concentrations.  It  Is  just  this  situation  that 
was  observed  In  Test  14A. 

No  data  was  available  for  cabin  materials  at  flux  levels  above 
4.5  Btu/ftz-sec.  The  current  state  of  knowledge  of  solids  combustion  does  not 


suggest  a  trustworthy  way  of  estimating  the  required  information.  Therefore  no 
attempt  was  iuads  to  improve  the  extrapolation  procedure.  Clearly,  more  work, 
both  theoretical  and  experimental ,  needs  to  be  conducted  on  this  problem. 
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